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PLASTIC DEFORMATIONS DUE TO DYNAMIC LOADING 
OF A BEAM WITH AN ATTACHED MASS! 


By T. J. MENTEL? 


ABSTRACT 


The influence of an attached mass on the large plastic deformations of a beam 
subjected to dynamic pressure loading is examined. The particular cases worked 
out in detail are those of a simply supported and a built-in beam, subjected to a 
uniform pressure pulse of a rectangular shape. In the analysis elastic strains are 
assumed to be negligible. The results contrast with those of the uniform beam 
in certain respects. The successive configurations of deformation excited by 
loads of various magnitudes are more complex than in the uniform beam (zero 
attached mass) case. Moreover, a new phenomenon appears, namely that plastic 
strains in the beam adjacent to the attached mass take place first in one direction 
and subsequently in the opposite sense, when the load magnitude exceeds a 
certain value. This reversal of sign of plastic flow may be important from the 
point of view of rupture, although the general effect of the attached mass is to 
strongly reduce the magnitude of the deflections of the beam. A discussion of 
the limitations and possible accuracy of an approximate solution based on a 
simple one-dimensional model is included. 


INTRODUCTION 


The primary purpose of this investigation is to study the effect of a mass 
attached to a beam (Fig. 1) on the deformations produced by dynamic pressure 
loading. The results are also applicable at once to the case of a portal frame 
(Fig. 2) in which the horizontal member is assumed to be rigid and to have 
negligible depth, and in which the load is assumed to be applied equally to 
both columns. 

The so-called plastic-rigid type of analysis, which is used in this paper, has 
already been applied to several problems of uniform beams involving various 
end conditions and types of dynamic loading (8, 4). The basic assumptions of 
this method, as applied to transverse deformations of beams, can be sum- 
marized as follows: 

(i) Only plastic flexural deformation is significant, and this occurs only 
when the moment at any cross-section reaches and maintains a fixed magnitude 
Mo, termed the “limit moment”’ of the beam section. 

(ii) ‘‘Plastic hinges” occur at all sections where Mp) is maintained. When the 
value M, is maintained at a given cross-section of the beam, the plastic hinge 


1Manuscript received February 1, 1955. 

Contribution from the Division of Applied Mathematics, Brown University, Providence, 
Rhode Island, U.S.A. The results presented in this paper were obtained in the course of research 
sponsored by Norfolk Naval Shipyard under Contract N189s—92568. 

*Research Assistant in the Graduate Division of Applied Mathematics, Brown University. 
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Fic. 1. Beams analyzed in paper. 
Fic. 2. Equivalent portal frame. 
Fic. 3. One-dimensional model. 


allows a discontinuity in the slope of the beam to occur across this section. 
If the moment has magnitude M, over a finite distance, the plastic ‘“‘hinge’’ 
becomes a finite plastic zone. 

(iii) Shear and axial forces occurring at plastic hinges are not large enough 
to have an appreciable effect on the value of Mo. 

(iv) All displacements are assumed small so that changes in beam geometry 
can be neglected in setting up the equations. 

This concept of localization of deformation at plastic hinges has led to useful 
methods for analyzing the failure of beams and frames under static loads (1, 5). 
A criterion has been offered (2, 6) for the estimation of the validity of the same 
concept in dynamic problems, where the neglect of elastic vibrations is in- 
volved. This criterion is based on the hypothesis that a plastic-rigid type of 
analysis will give good results (as far as the neglect of elastic strains is con- 
cerned) when the work done in the plastic deformation greatly exceeds the 
maximum elastic strain energy that can be stored in the beam. This can be 
expressed as a requirement that the load (expressed nondimensionally) should 
sufficiently exceed a certain number which is a constant times the square of the 
ratio of the duration time of the load to the fundamental period of the beam. 
However, for this present analysis, we will accept the assumptions of the 
plastic-rigid theory and not discuss their validity further. 

A device frequently employed in studies of dynamic loads on structures is 
the representation of these structures by semirigid models which consist of a 
finite number of rigid masses interconnected by springs with proper character- 
istics. The simplest case is of course the one-dimensional model (Fig. 3) which 
might be taken as a working approximation of the portal frame of Fig. 2. 
This approximation might be expected to be valid for certain of the extreme 
cases, such as where the attached mass is very small, or where it is very large 
in comparison to the remaining mass. However, in the intermediate cases, it 
is not obvious how the masses should be lumped in order to obtain reasonably 
correct results, if indeed reasonable approximations are feasible. One purpose 
of this report is to give information concerning the accuracy to be expected 
from solutions obtained by using such a one degree of freedom model. 
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FORMULATION OF GENERAL DYNAMIC EQUATIONS 


We consider the effect of a uniformly distributed load W(t) on a simply 
supported beam of constant limit moment Mo, constant mass m per unit 
length, and carrying a concentrated mass Q at its center (Fig. 1). Since it is 
found that only slight changes occur in the formulae in treating the case of a 
built-in beam, only the case of simple supports will be discussed in detail, the 
results for the other case being quoted at the end. 

In formulating the kinematical behavior of the beam under the action of 
the load, we recall that in the plastic-rigid analysis no deformation takes place 
until the load is such that plastic action can be sustained. We have, therefore, 
an initial period during which the load builds up (if the load is assumed to start 
from zero or some small finite value) and the beam is assumed to remain rigid. 
The problem at this stage is statical, and we readily solve for the maximum 
moment in the beam, which in this case occurs at the center and attains the 
value My when W = 4M,/I. The first phase of motion to occur, then, as sug- 
gested by the statical moment distribution, would seem to have a plastic hinge 
at the center (or two hinges infinitesimally close together on either side of the 
mass). Upon postulating such a condition, however, it is noted that for subse- 
quent deformation to take place, a shear force would be required at the 
concentrated mass in order to bring it into motion. Such a force cannot be 
allowed, for, as shown in Fig. 4, it would demand the occurrence of a moment 
greater than M, at some point within the beam. It is therefore concluded that 
if a hinge appears at the center, it must immediately travel outwards, and for 
this reason, we consider the first phase of motion to be characterized by two 
moving hinges as shown in Fig. 5. 


Vertical 
shear 
force 
diagram 





F [Mo 
Bending b 
moment 
diagram Fig. 4 


8 Plastic hinge 





Fic. 4. Formulation of first phase. 
Fic. 5. Deformation mechanism for first phase. 
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The resulting motion can then be described by the following dynamical 
equations: 


linear impulse-momentum, 
t t ail dy ds 
[1] wat-2 | Rat = 2 f m\ — } dx+{2m/[1—&]+Q} — ; 
J t, i 0 dt dt 
vertical acceleration of outer portion of — 


[2] Wt-2R = me ee — 


angular acceleration of outer portion of oe 


Wel yy, = met dO 
[3] 7 Me = a 
vertical acceleration of center portion of beam, 
, d’s 
[4] W(1—§) = {2mi[1—§]+ Qi} 7a ; 


where (see Fig. 5) 


Ly cartesian coordinates, 


> 


rotation of outer beam segment, 


wr 


— dimensionless parameter, 
! — beam semispan, 


~ 


transverse displacement of attached mass, 
R — reaction at supports, 
to — time of start of beam deformation. 


Now defining » = W//Mo, K = Q/2ml, n = t/r, where 7 is a characteristic 
time of the force pulse, and letting the dot signify differentiation with respect 
to n, we can reduce the above equations to 

“af 3\ ml” 

J yutt— 1-26 dy = 3 ajo et 2K}, 
e 


No MoT 


(6] jl (1—&+K) = 5(1 - 


mor 
In evaluating the integral on the right-hand side of equation [1], we have used 
the following condition of continuity of transverse velocity: 
do tl ds 
dt dt ° 


Finally, a single equation is obtained by differentiating equation [5] and 


solving with equation [6] 


7] G8 2K) al ~ 4] hee Pfau a3} dn 


2_u(1—£) _ 
2(1-E+K) ° 


ic 
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A different form of this equation can be obtained by using the jump in trans- 
verse acceleration across the traveling hinges as developed in (5). In this case 
we obtain 
= 6&l + 6é/ | 
where 5 and 6 are given by equations [3] and [4]. The resulting equation, after 
elimination of 6 by differentiation, is 


| w(l-t) 3 a wk 
p) (AG aO. 3 oat afer eS oto 
Bhat us @ 7 
Haas taK ie 0. 


[t appears that, in general, some iterative numerical procedure would have to 
be used in order to solve either of these equations. However, in considering 
forms of the pulse u(/), we note that previous work on similar problems (3, 4) 
has shown the total deformation to be largely insensitive to the type of pulse, 
and hence the selection of a rectangular pulse, which offers considerable 
analytic simplification, weuld seem to be a reasonable practical choice. In 
this case we have 
a = & = O during the pulse (0 <t <1), 
uw = O after the pulse (7 < ¢) 


Equation [7] can then be rewritten 


(9 ] Ho — 120 —f+K) 

~ &>) (1—f+3K) * 
Equation [9] can be used to determine the initial position of the hinges for 
pulses starting with a finite value » > 4, as well as giving the position of the 
hinges during the action of a rectangular pulse. Fig. 6 shows a plot of this 
equation for the full range of the parameter K. 

For K approaching zero, which represents the case treated in (4), we note 
that the hinges remain at the center of the beam for values of the load 
4 <u < 12and spread into finite plastic regions for higher loads, in accordance 
with the results obtained therein. We also note, in the present problem, the 
appearance of a new phase for higher loads, where a negative plastic moment 
is developed at the center because of the inertial resistance of the concentrated 
mass. The criterion for the appearance of this next phase can be obtained by 
solving for the moment at the center of the inner beam section, which is 
illustrated in Fig. 7. Setting the moment at the center equal to — Mo, we obtain 


(10] §=1- 
The intersection of the curves defined by equation [9] with those of equation 


[10] is shown by the broken line in Fig. 6, which thus defines the upper bound- 
ary of the first phase. 
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Fic. 6. First phase hinge locations. 





Ler . 


Fig. 8 


Fic. 7. Inner beam segment as a free body. 
Fic. 8. Deformation mechanism for second phase. 


The mechanism of the second phase is shown in Fig. 8. Using the same 
procedure as before, we obtain the following equations of motion: 
linear impulse-momentum, 


3 Soe 
(11) f{ou-1-3 dy = fs O12] tall —2+E]4+=F° 5 
continuity of transverse velocity, 
[12] § = 6&1 + a(l — &)/; 
jump in transverse acceleration, 
[13] S = bf] + £1(6 — &) + a(l — é)/; 


acceleration of concentrated mass, 





ume aan I AI 


ne 
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0 _ mt  __ Mta(i—s) 
| “Mor” = -2(1—#)(1—£+4K)’ 


angular acceleration of inner beam segment, 


ml’ _ 24(1—§-+K)—3yK (1—8)" 


[15] “Me (1—¢)°(1—£+4K) ; 
and angular acceleration of outer beam section, 
. mi” 3 2 
6 —— a 7 "—4), 
(1 i] 6 Mor’ 4¢° (ué ) 


We can now formulate the general equation of motion by using either 
equation [11] or [13] as a basis, and employing equations [12], [14], [15], and [16] 
as required. Either of these formulations, however, results in an unwieldy 
non-linear equation which does not appear to lend itself to analytical treat- 
ment. For the special case of the rectangular load pulse, we can use equations 
{13}, [14], [15], and [16] to obtain 
nt _ (12/0") (1 —fo)"(1— fo + 4K) — 48(1 — fo + 2K) 
| 7] Ho = 3 : 

(1—£) 
A plot of the curves defined by equation [17] is shown in Fig. 9. It is noted that 
during this phase ) becomes independent of Mp, for infinite A, which corre- 
sponds to the case of a uniform beam simply supported at one end and built 
in at the other. 





HA tant KI} 
= 
—K#10}& 
K+ 
Plastic regions 
of third phase 


Upper boundary 
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Fic. 9. Second and third phase “hinge” locations. 
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A third phase can be identified by noting the conditions under which the 
plasticity condition is violated in the second phase, i.e., when a moment 
greater than Mp, is developed at some point in the beam. If we consider the 
outer beam segment first, we see that the plasticity condition might be violated 
only at the hinge end, since the moment must be zero at the simple support 
and increase on going inward. The condition that the moment decreases 
steadily from M, on going outward from the hinge is obtained from the require- 
ment that d?M/dz* is of one sign in this segment. The following condition 


ensures this: 
[18] mt6 < Wr?/2/. 


The critical value is found by using the equality sign in equation [18], and 

this with equation [16] reduces to 

[19] he? = 12, 

where the subscript is now used to denote the particular plastic-rigid interface. 
In the case of the inner beam segment, we have a moment distribution as 

shown in Fig. 10. From the shape of this moment curve it is clear that no 

violation of the plasticity condition can occur near the inner end of this seg- 

ment. Such a violation would, however, occur at the outer (traveling hinge) 

end of the segment if d*\/7/dz? changed sign near this end. The condition that 

no violation occurs is 


[20] may < Wr?/2l, 
ao being the acceleration of the outer tip of this section. Now we have 
[21] a = 5 —a(l — &)/. 


Using equation [20] with an equality sign, together with equations [14], [15], 
and [21], we obtain the criterion 


[22] w= 24 ( 1—&+K) —12(1-&) 
oe K(1 $s y e 


A plot of the curves defined by equation [22] has been included in Fig. 9, 
where their intersection at the same point with the curve defined by equation 
{19] and the appropriate curves of equation [17] furnishes a check of our 
assumed deformation mechanism. 

The third, and now necessarily final phase, thus has two finite plastic zones 
as shown in Fig. 11. An initially infinite load thus would imply that the entire 
beam is instantaneously plastic. The boundaries of the plastic zone are defined 
by equations [19] and [22] for the case of the rectangular pulse during the 
time of its action or for an increasing load. After such a pulse, or for a de- 
creasing pulse, equations [19] and [22] become only inequalities and either 
impulse-momentum equations, such as equation [11], or an acceleration jump 
equation, such as equation [13], are required to describe, the motion. However, 
we shall not carry through the analytical treatment for this phase, and will 
limit ourselves to maximum values of uo below about 70 so that only the first 
and second phases occur. 





id 
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Fic. 10. Moment distribution for inner beam segment. 
Fic. 11. Deformation mechanism for third phase. 
Fic. 12. Deformation mechanism for final motion. 


BEAM DEFORMATION 


The deformation of the beam is best described by the transverse displace- 
ment of the concentrated mass and the angular displacements at the ends and 
at the center. Considering the transverse displacement first, the complete 
displacement consists of the following parts: 

(1) Displacement s; during the time of action (m: — mo) of the load pulse. 

(2) Displacement s2 during the interval (n2 — m) when the motion corre- 
sponds to the kinematical picture of Fig. 11. This interval ends when the finite 
plastic zones are reduced to point hinges. 

(3) Displacement s; during the interval (n3; — n2) when the motion corre- 
sponds to the kinematical picture of Fig. 8. This interval ends when the 
angular velocities of the inner segments vanish. 

(+) Displacement s, during the interval (ns — 3) when the motion corre- 
sponds to the kinematical picture of Fig. 5. This interval ends when the 
traveling outer hinges reach the center of the beam. 

(5) Displacement ss; during the interval (ny; — 74) when the motion corre- 

sponds to the kinematical picture of Fig. 12. This interval ends when the 
angular velocities of the half-beams vanish. 
For the general case of a load starting from zero and gradually increasing, 
the displacement s; is subdivided into three displacements occurring during 
the three phases of motion. In the case of loads initiating only the first and 
second phases, s2 is zero, and when only the first phase is initiated, sz and 83 
are zero. 

For rectangular load impulses (which we are considering) the displacement 
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s; is obtained by direct integration of equation [6] for the first phase, and of 
equation [14] for the second phase. The resulting expressions are: 
load producing first phase 

ml? __wo(1—Eo)_ 

“' Mor 4(1—fo+K)’ 


load producing second phase 


ml* a 24+ (1 — £0) 


* Mor’ ~ 4(1—&)(1—&)+4K) 

Assuming the load magnitude to be such that s. = 0, we have next to find 
the displacement s; by solving either equation [11] or [13] with » = 0 and ¢ 
being variable. Because of the non-linearity of these equations a numerical 
integration would be required to obtain exact results. However, an approximate 
solution can be obtained by recalling that during the second phase ~ remained 
essentially constant for the smaller values of load. We therefore inquire into 
the value of ~ during the second interval 7; — m to see if a similar behavior 
occurs. We will restrict ourselves to values of uo not greater than 50, and asa 
first approximation, will calculate the value of & at the beginning and end of 
the interval for the same value of &. 

Using equations [13], [14], [15], and [16] we obtain 


ml _ 12¢°(1—§+2K) —3(1—£)*(1—£+4K) 


23] g tio ia ese. es 
Mor &(1—£&)"(1—E+44K) (0—a@) 
where at the beginning of the motion 
ml" 3 __24(1—fo+K) —3uoK (1 — fo) 


/ Be eee 
¥ (Moko +) (1 —£)*(1 —f)+4K) 


The approximate value of the angular velocity difference at the end of the 
motion is 


(0-4) an = Ge, 


ml _ 3u0S, 2 _K(1—&)" \ 


(0—a) Mor 4te™? ~ 2(1—&)+-K)) 


where we have used the value 


woK (1—£0)” 


» : = 
[24] 1 8(1—f+K) 


which is obtained by setting & = 0. 

For large A, the numerator in equation [23] tends to zero while the denomin- 
ator remains finite; thus for A = 10, wo = 50, the approximate initial value of 
— is 0.003. For K = 1, wo = 50, the approximate initial and final values of 
€ are 0.03 and 0.07, with the approximate value of (n3 — q2) = (nx — m) being 
0.4. Only a small variation in & is therefore indicated, and comparison with 
results obtained by an iterative numerical procedure using equations [11], [14], 
[15], and [16] shows that for this range of yo, the assumption of constant £ 
during the second interval gives final deformation values within three per cent 
of the exact values. 
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Using this approximation of constant £, we readily obtain the following 











relations: 
(25] ¢, TE = — 24m uo(1— fo)" _ 
, “* Mor’ 2(1—£&)(1—0+4K) ’ 
(26) _, mm _ 12 (ms'=1) +mo(1 = f0)*(9s— 1) 
o "* Mer’ 2(1—&) (1—f0+4K) : 


where 7; is defined by equation [24]. The velocity §;, of the concentrated mass 
at the end of the second motion, has been included because of its usefulness 
in computing the remaining displacements. 

For the motion in the next interval, (7, — 7:), we set u» = 0 in equation [5] 
to obtain 


s1( me Pr ” 
ef sotto 2y dn— [Ban = 5 (2—£42K). 


0 


We proceed by differentiating equation [27], noting that equation [4] shows § 
to be constant during this type of motion. Subsequent integration then yields 


ef ml _ 3 
I Mr ht? 


where c is evaluated at 7 = 1, & = &. The time 4, determined by setting 


§ = |, is then 
ie , mt. it 
i = nat+53 uel 6 . 


The corresponding displacement can then be written as follows: 


ml \; epee) 
“Mer” Med V6 OS: 


The mechanism of the final motion is illustrated in Fig. 12. Moment equi- 
librium of the half-beam gives 





- ml —.. a 
“Mor (1+3K)° 


Integrating once and evaluating the constant at » = m, § = $3, we obtain 


. a ml’ (143K 
n—m = (8:—3) Me\ 3 )° 


The motion ceases when § = 0, or 


. ml (143K 
s— 14 = S3 Mor —— 


Performing a further integration, and using the conditions that when » = m, 
s = 0 and when 7 = 75, s = 55, we get 


ml” f ml” ee 


8M = 6 


[28] 


"8 Mer 
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which can be reduced to 


i molt) Vf 143K} 
“Me (G—-&+K)) ( m4 )° 


The total transverse displacement can then be written 





[29] , nl _ 24ns°+uo(1—£o)*(2n3—1) i aff eaten 
oS Mor ~~ 40) fe 4K) 8 Mor? 6 


where n; and §3 are given by equations [24] and [25]. For the case of load pulses 
initiating only the first phase, the expression reduces to 


ml? _ _wo(I—fo)_J1, _wo(1—E0) ae 
2(1—&+K) 6 


A plot of equations [29] and [30] is shown in Fig. 13. In the special case where 
K = 0 and & = 1 simultaneously, equation [30] can be evaluated only by 
using a limiting process. A simpler procedure is to note that the proper ex- 
pression to be used for the velocity § is easily obtained from equation [16] by 
putting = £ = constant; the result is 


BO) sa = 2(1—&+K) \2 


en Deta*.. 
Ss Mor’ fer Ato” (uoko 4). 
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Fic. 13. Displacement at center due to a uniformly distributed rectangular load pulse 
(simply supported beam). 

Fic. 14. Angular displacement at end due to a uniformly distributed rectangular load 
pulse (simply supported beam). 
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Using this expression with & = 1, equation [30] becomes 


ml 3 
[31] ‘Ma * 35 Ho(mo—4), 


as obtained in (4) for the condition of K = 0 and 4 < wo < 12. 
The angular displacement at the ends of the beam is determined in a manner 
completely analogous to that of the transverse displacement, the result being 


o ott « Bias A | aay 126+ 3K) 
[32] Ot = amt {2n3—1]—4n3 +53 Mor2h \ 6 


When the second phase does not occur the expression is: 


a 3 _Ho(1— £o) ep Vian 


A plot of equations [32] and [33] is given in Fig. 14. 

For the case of K approaching infinity, the behavior tends to resemble that 
of two separate uniform beams, and thus gives an insight into the formation 
of the plastic-rigid deformation mechanism for a uniform beam built-in at one 
end and simply supported at the other, when the built-in support is allowed to 
be displaced. For infinite K, the resemblance to the case of two uniform beams 
is exact, and no deformation occurs until the load reaches the critical value 
which initiates the second phase. This value of the load is easily found to be 
uo = 23.315 with the corresponding value of & being 0.4142. 

In addition to the transverse displacement at the center and the angular 
displacement at the ends, an important measure of damage for certain ranges 
of uo and K may be the angular displacement at the concentrated mass. It is 
interesting to note that when the load is high enough to produce the second 
phase, yielding occurs under negative moment (0 to 73) followed by yielding 
under positive moment (4 to ns). This reversal of the plastic flow in a problem 
of dynamic loading of a structure appears to have been observed here for the 
first time. It may have important significance in connection with the possibility 
of rupture of the beam. 

The angular displacement at the concentrated mass is described by the 
following expressions: 


(a) for loads producing the first phase 


mi _f§ (1—&) \*f eae 
“Mr ~ \A—-&+K)) \ i 

and the time variation is given by 

ml’ _ _(1—f) theca wail 

® Mor ~ 21—&+K) "~™ —20+3K)"—™ } 
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(6) for loads producing the second phase, 
ml _ 24(1—f-+K)—3woK (1 —£0): 














“Me” 2i-—byd-btn) 
a, = 0, 
‘ ml" oe —3{uoK (1— 0)" —8(1—fo+K)}? 
*Mor  16(1—£0)*(1—f0 + 4K) (1—f0+K) * 
a, = 0, 

mE { pelt) ieee 
“8 Mer U—te+K)) V 24 


A plot of the above results is shown for two cases in Fig. 15. In the first case, 
the load is just under that for the second phase, while in the second case, the 
second phase is excited and the reversal of plastic flow is displayed. 
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Fic. 15. Time variation of angular displacement at the mass. 


The deformation mechanism for the first phase in the case of a built-in 
beam is shown in Fig. 16. The only changes which occur in the dynamical 
equations are the values of certain constants. For this reason, all intermediate 
steps in their development are omitted, and only the final equations are shown 
below. 
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ria) f{utt—ae1—8} dn = dee epat—2er 4 2Rsh 





_ (24/ée° ) (1-0)? (1 = fo 4K) — —48(1— fot 2K) 


[174] if)" 
‘ a _ml" - Ses: 
284 ] Me Te’ 


The corresponding expressions for the total deformations are 


mi? 24ns?-+-uo(1—Eo)?(2ns—1) + ;, mt yaa 
"Me |0=—- 40-8) 4a) * Mor’ 12 , 


204] sy ee 


[29A ] 











Mor” a 2(1—f)+K) 12 
i aml *59—28,4-8K 





saa 0! yo(1—£o) \"fa-2 mte8) 


2 
9 or? = gf (uote — “evefgmt 2(1—%+K)S 12 


The values of 73 and §; are given by equations [24] and [25] respectively. A plot 
of equations [32A] and [33A] is shown in Fig. 18. 
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Fig. 16 Fig. 17 


Fic. 16. Deformation mechanism for the first phase of the built-in beam. 
Fic. 17. Deformation mechanism for the first phase of the uniform beam. 


USE OF A SIMPLE MODEL 


A standard technique employed in engineering analysis of structures under 
dynamic loading is the replacement of the structure by a model having a finite 
number of degrees of freedom. The simplest such model is the single mass- 
spring system (Fig. 3) of one degree of freedom. As noted in (4), if in our 
problem we assume the beam to be displaced in a definite deflection curve, 
we reduce its infinite number of degrees of freedom to only one, thus making 
a single degree of freedom model the appropriate simple counterpart. Under 
these conditions, the properties of the model can be chosen so that the dis- 
placement of the model corresponds exactly to the deflection of, for example, 
the midpoint of the beam. However, if (as in the beam problem considered 
here) the deflection curve changes its shape during the course of the motion, 
the one-dimensional model can give exact results only during that part of the 
motion in which the selected mode is maintained. 

One procedure for reducing the real structure to the simple model, once the 
form of the model has been decided upon, is to relate the applied loads, which 
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in our case are the distributed load for the beam problem and the concentrated 
load for the model, so that the work done is the same in each case; and to 
choose a mass for the model so that its total kinetic energy is the same as that 
of the beam moving in the assumed deflection curve. 

In the case of a simply supported uniform beam, for values of the load in 
the range 4 < wo < 12, a plastic hinge occurs only at the center of the beam 
(Fig. 17). There is therefore only one mode, and we obtain exact results by 
using a simple model. Following the procedure outlined above, we equate the 


- - - Approximate solution / 


—— Exact solution 


First phase 
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Fic. 18. Angular displacement at end due to a uniformly distributed rectangular load pulse 
(built-in beam). 

Fic. 19. Comparison of approximate and exact solutions for the displacement at the center 
due to a rectangular load pulse (simply supported beam). 


work done by P and W to obtain P = W/2, and similarly equating the kinetic 
energies we find M = 3ml. Lastly, we equate the works done by the resistive 
forces to obtain F = 2M,//. Now for the simple model, where we allow the 
spring to offer plastic resistance only, the expression for the total displacement 
resulting from a rectangular load impulse of magnitude P and duration 7 is 
34] M_Pijp_,\ 
3: Say et See he 

Fr 2F\F US 
Substituting the relationships determined for P, M, and F into this equation, 
we obtain 


. ml” 3 
[35] S Mort = 39 Ho(uo— 1) 


which is identical with the exact expression (equation [31]) as furnished by a 
direct solution of the beam problem for yo in the stated range. 
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From the complete solution of the beam problem (4), we know that for 
uo > 12, the appearance of a finite plastic zone introduces a new kinematical 
behavior, which is reflected in the expression for the deformation. However, 
had we proceeded only with the model, we would have had no means of 
accounting for this phase, and Fig. 19 has been included to show the error 
incurred by using the elementary expression beyond its range of validity. 

Considering now the problem of a beam carrying a concentrated mass, one 
might expect that the replacement of the actual system by a one-dimensional 
model would be more satisfactory than in the case of the uniform beam. The 
natural beam configuration to assume in the reduction to the one-dimensional 
system is again the configuration of Fig. 17, except that the central “‘hinge”’ 
now consists of two hinges adjacent to mass Q, as in Fig. 12. Using the same 
procedure as before in obtaining the correspondence between variables, the 
only change found to occur is in the expression for the mass M of the model, 
which now becomes 


M = 3ml+ Q. 
Substitution into equation [33] yields the result 
}? 3 
9 ey i aa yy). 


Mor ~—-32(1+3K) 


A comparison of this result with the exact results for two values of K is also 
shown in Fig. 19. The error is noted to be relatively small and positive for 
K = 0.2, and, contrary to expectation, to become increasingly large and 
negative for large values of K. 

In an attempt to improve the approximate results, a second choice for the 
configuration might be that shown in Fig. 5. With this configuration, the 
equivalence procedure with the one-dimensional model leads to the following 
relations: 

P = WIR — &)}, 
M = ml{3(6 — 4&)} + Q, 
F = 2M /éd. 


Using equation [34], the new expression for the total deformation is 


J ml _ __3(2—0) mo 
Mor’ = 32(3—2&)+3K) 





(uol2—£o] o—4). 


However, a numerical analysis shows that, although some improvement can be 
effected by choosing £) to be approximately 0.9, it is negligible in comparison 
to the total error. 

We have already noted that as K becomes large, the behavior during the 
deformation tends to resemble that of two uniform beams built in at the center. 
This might suggest the alternative approach of actually treating the problem 
as two uniform beams for the case of large K. However, since the deformation 
obtained by these results would not commence until uo attained a value of 
23.315 for the simply supported beam, and 32.0 for the built-in beam, it could 
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hardly be accurate for values of wo of this order for K ranging from say | to 5, 
for at these values of wo and K, Figs. 13, 14, and 18 already show appreciable 
deformation. We note, of course, that for both large u» and K, the most im- 
portant parameter is the angular deflection at the ends. However, the model 
we have used is not directly amenable to duplicate angular displacement, and 
the simple expedient of dividing its answers by / leads to grossly inaccurate 
results for large K. We might, alternatively, consider a different model which 
incorporates angular instead of rectilinear motion, but it is evident that for 
values of K ranging from 1 to 5, and wo of the order 30, the generally radical 
changes experienced by the deflection curve could not be effectively duplicated 
by any simple one-dimensional model. A better account might be made of the 
original continuum problem if further characteristics were added to the model, 
but this would seem to defeat the essential purpose of simplicity for which 
the model was introduced. 


CONCLUSIONS 


A simply supported or built-in beam, carrying a concentrated mass at its 
center, may experience as many as seven distinct patterns of deformation 
during its whole deformation process resulting from a uniformly distributed 
load pulse. The first effect of a small but finite mass at the center is to prevent 
the development of a centered plastic zone which appears in the second de- 
formation phase of a uniform beam. As the relative size of this mass is in- 
creased, a new second phase becomes important. This second phase is char- 
acterized by an additional plastic hinge at the center where the direction of 
deformation, as shown in Fig. 8, is opposite to that which occurs in the uniform 
beam case. When the ratio of the attached mass to the beam mass becomes of 
the order of ten or greater, the deformation behavior closely resembles that of 
two uniform beams, with their inner supports being built-in. 

A second effect of the attached mass is the occurrence of a reversal in the 
plastic flow at the middle of the beam when the load is high enough to initiate 
the second phase. This behavior may have an important influence on the 
tendency of the beam to fracture and hence the present results may provide 
a useful parameter in formulating an experimental criterion for beam failure. 

The technique of approximating a continuous structure by a simple model 
in dynamical problems provides good results only when the motions to be 
accounted for are nearly identical, so that they may be closely described by 
using only a small number of parameters as provided by the simple model. 
In the case of a uniform beam acted upon by a small transverse load, only one 
deformation mode occurs, and hence a one-dimensional model can be con- 
structed to duplicate exactly any chosen characteristic of the deformation. 
However, if a concentrated mass is attached to the beam, or if high values of 
load are employed, marked changes may occur in the mode of deformation 
and a simple one-dimensional model in general gives erroneous results, though 
fair approximations may possibly be obtained for restricted ranges of the 
parameters involved. 
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_SOME CHARACTERISTICS OF SPRAY-DRIED WHEAT GLUTEN! 


By W. B. McCoNNELL 


ABSTRACT 


Dispersions of wheat gluten in 0.01 N acetic acid were spray dried using 
various drying conditions to give powders retaining most of the original baking 
strength of the gluten. The dispersions had a pH near 5.4 and for most experi- 
ments they contained from 8 to 9% solids. With the spray drying equipment used 
an inlet drying air temperature of 160°C. and outlet temperature of 90°C. 
consistently gave gluten which compared favorably with that obtained by freeze 
drying. Baking (100-gm. loaves) data obtained after addition of spray-dried 
gluten to a cake flour gave a regression coefficient for loaf volume on per cent 
protein added of 48.2. The preparations contained 15.4+0.3% nitrogen and 
could be readily redispersed in 0.1 N acetic acid. The moisture content of the 
freshly dried gluten preparations ranged from 0.4 to 4%; however they adsorbed 
atmospheric moisture and at 90% relative humidities they became wet and 
gummy. The preparations were uniform, consisting largely of spherical powder 
particles with diameters between 10 and 20 microns. It is concluded that it is 
technically feasible to prepare and use spray-dried gluten for improving the 
baking quality of low protein flour. 


INTRODUCTION 


In 1950 Lusena (3) described a method for preparing a dispersion of pure 
wheat gluten in dilute acetic acid solutions and showed that the gluten ob- 
tained by freeze drying the dispersions was essentially undenatured. Pre- 
liminary work by Lusena and Adams (4) indicated that the gluten dispersions 
could be spray dried to give a finely powdered product retaining much of its 
bread-making quality. This work has been extended to evaluate the practica- 
bility of spray drying gluten for use as a flour additive in the baking industry. 

This communication reports on the nature and baking characteristics of 
spray-dried wheat gluten. It includes some observations on the effects of 
variations in the drying process on the quality of the product and gives data 
related to physical characteristics of significance when considering its possible 
production on a commercial scale. 


EXPERIMENTAL 
Preparation of Gluten Dispersions 


Acetic acid dispersions of gluten from an unimproved top patent flour 
(12.8% protein) were made according to the method outlined by Lusena (3). 
A laboratory homogenizer (Eppenback Homomixer No. 2) was used to dis- 
perse the gluten in 0.005 N acetic acid. Four minutes’ agitation at full speed 
was sufficient for thorough dispersion of the gluten from 1500 gm. flour. The 
gluten was precipitated with calcium hydroxide, and redispersed for drying 
in 0.01 N acetic acid (3). The second dispersion had a pH near 5.4 and con- 
tained from 8 to 9% solids. The nitrogen content of the dry solids was 15.4%. 


1Manuscript received February 4, 1956. 

Contribution from the National Research Council of Canada, Prairie Regional Laboratory, 
Saskatoon, Saskatchewan. Issued as Paper No. 189 on the Uses of Plant Products and as N.R.C. 
No. 3622. Presented in part at the Annual Conference of the Chemical Institute of Canada at 
Windsor, Ontario, June 1958. 
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Spray Drying Method 

The gluten dispersions were dried in an insulated stainless steel spray drier 
with a maximum inner diameter of 24 in. It was constructed according to the 
detail of Woodcock and Tessier (7). 

In operation, the drier was first fed with acetic acid until temperature 
equilibrium at the desired conditions was attained. After the gluten dispersion 
was admitted through the spray nozzle minor adjustments were usually neces- 
sary to maintain the initial temperature conditions. After drying was com- 
pleted the feed lines were flushed again with dilute acetic acid. 

A number of preliminary experiments were necessary to obtain a satis- 
factory adjustment of the atomizing nozzle and atomizing air pressure. If 
the spray was coarse some undried droplets reached the inner surface of the 
cone where they were deposited as a hard coating, adhering firmly to it. The 
same difficulty was encountered with a wide-angle spray. Wide sprays also 
increased deposition of gluten on the spraying nozzle with a resultant dis- 
tortion of the spray pattern and eventual clogging of the nozzle. Best results 
were obtained from a fine spray with the droplets directed well toward the 
apex of the cone. An atomizing air pressure of 25 p.s.i. was used in all reported 
experiments. 

About half of the dried gluten collected as a fine white powder in the re- 
ceiving jar, and in the bottom of the cone. This was designated as Fraction A. 
A considerable amount of material adhered loosely to the upper portion of 
the walls and was collected by brushing and rubbing the lining of the drier 
(Fraction B). 

Comparison of the two fractions was made throughout the whole series of 
experiments. In general Fraction A was only slightly superior to Fraction B 
in baking quality, but it was considerably more uniform in particle size and 
had a greater bulk density. The material which adhered to the walls of the 
drier came, at least in part, from droplets that reached the sides of the cone 
before they were completely dry. Since this difficulty would likely not be 
encountered with a larger spray drier the free powder was considered to be 
more characteristic of a spray-dried gluten. Results reported herein are there- 
fore limited to those obtained with Fraction A. 


Analytical Methods and Baking Tests 

Nitrogen determinations were done by the conventional Kjeldahl methods. 
Moisture contents of the preparations were obtained from the loss in weight 
of samples heated at 110° C. for 20 hr., and total solids in the dispersions were 
measured by evaporating aliquots to dryness. 

Samples of freeze-dried gluten dispersions were prepared with a laboratory- 
scale glass lyophil apparatus (6). 

The gluten quality was evaluated by observing its ability to improve the 
bread-making quality of a low protein flour. For the experiments in which the 
drying conditions were varied, a 9.3% protein pastry flour was raised to 15% 
protein by addition of the spray-dried preparations. The somewhat larger than 
normal final protein content was advantageous for detecting differences in 
gluten quality. 
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The gluten and the base flour were mixed dry and the resulting flour was 
test-baked, using a sugar-rich baking formula which maintains adequate gas 
production throughout the baking procedure* and thus permits the evaluation 
of flour strength. The volumes of the loaves were measured. The loaves then 
were scored for external characteristics and photographed to provide a perma- 
nent record. 


EXPERIMENTAL RESULTS 
Effect of Drying Temperatures on Gluten 


Table I contains a summary of data obtained to characterize gluten samples 
prepared by operating the spray drier with a series of different inlet and outlet 


TABLE I 


BAKING CHARACTERISTICS AND MOISTURE CONTENTS OF GLUTENS SPRAY 
DRIED WITH DIFFERENT INLET AND OUTLET AIR TEMPERATURES 





Fraction A 











Outlet Inlet % N in Av. loaf volume, 

Exp. temp., temp., dry Moisture, cc. (pastry 

No. PS. ° solids % flour base) 
1 80 170 15.3 1.85 795 
2 86 160 15.5 1.20 750 
4 80 140 15.6 1.68 750 
5 80 130 15.1 2.70 730 
6 90 170 15.5 3.66 835 
7 90 160 15.2 1.22 845 
8 90 150 15.3 2.23 715 
9 90 140 15.7 0.78 735 
10 100 170 15.2 1.58 710 
11 100 160 15.2 1.14 740 
12 100 150 15.6 2.65 775 
13 100 140 15.6 0.83 740 
14 110 170 15.5 0.40 725 
15 110 160 15.4 1.35 770 
16 110 150 15.3 Lid 720 
17 120 180 15.1 1.36 540 


The average loaf volume given by the pastry flour (9.3% protein) was 496 cc. 
A standard household flour (12.8% protein) gave an average volume of 760 cc. 

An average loaf volume of 785 cc. was obtained by addition of freeze-dried 
gluten in the same way as the spray-dried gluten was added. 


drying air temperatures. Although there is a good deal of variability in the 
moisture content of the products, the values are in general agreement with the 
value of about 2.5% reported by Lusena and Adams (4). Most of the prep- 
arations contained less than 2% moisture and only one sample exceeded 3%. 
No consistent relation existed between the moisture content of the products 
and any of the variables measured. 

Throughout the whole series of experiments on A and B fractions the average 
difference between the volumes of duplicate loaves was 38 cc., and since most 
of the volumes reported in Table I fall within a 100 cc. range, the method as 


*One hundred grams flour, 4 gm. skim milk powder, 3 gm. fat, 3 gm. yeast, 5 gm. sugar, 1.75 gm. 
salt, 0.1 gm. NH,H2PO,, 1 ml. malt solution (60 gm. non-diastatic mali in 200 ml.). 
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used is evidently limited as a means for distinguishing between the prep- 
arations. 

It is apparent however from the data given that gluten retaining an apprecia- 
ble amount of its baking quality was obtained over a fairly wide range of 
conditions. A product comparable to freeze-dried gluten was obtained by 
operating the drier with an inlet temperature of from 160° to 170° C. and an 
outlet temperature from 80° to 90° C. Conditions within this range consistently 
gave a good product and were used for 12 preparations giving loaf volumes in 
excess of 800 cc. The most severely damaged sample. was made with an inlet 
temperature of 180°C. and an outlet temperature of 120°C. The trends 
suggest however that an inlet temperature of 180° C. could be used with an 
outlet temperature below 100° C. Results with these conditions are not re- 
ported because with the high flow rate required, severe clogging of the nozzle 
was encountered. It is significant, however, from a practical point of view, that 
some of the best products were made with the drier operating near maximum 
capacity. The deleterious effect of heat above 100°C. was also shown by 
heating freeze-dried gluten for one hour. Baking tests showed incipient damage 
at 100° C. followed by a linear decrease in loaf volume from 100° C. to 160° C. 
At the latter temperature destruction of baking strength was complete. 

Fig. 1 contains photographs of the loaves obtained from preparation 6 





Fic. 1. Comparison of loaves baked from household flour, pastry flour, and pastry flour 
with spray-dried gluten added. 
H—household flour; P—pastry flour; 6—pastry flour plus spray-dried gluten preparation 6. 
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(Table I) together with those prepared from household flour and pastry flour. 
Although the loaves shown were somewhat larger than average, they were 
typical in appearance to those baked with the added gluten and could not be 
distinguished in appearance from those baked with added freeze-dried gluten. 
The crust color and shape with preparation 6 was as good as with household 
flour. In general, loaf characteristics such as crust color and shape were some- 
what better for loaves of larger volume. All loaves containing spray-dried 
gluten were however slightly ‘“‘overaged’’. It should be noted that no oxidizing 
or reducing agents were added to the experimental doughs. 

Since the costs of preparing spray-dried gluten would depend on the concen- 
tration of solids in the dispersions, some experiments were made with different 
concentrations of gluten. Standard methods of preparing the solutions were 
used except that the final volume of dispersing liquid was varied to give the 
concentration desired. Drying was done with an inlet temperature of 160° C. 
and an outlet temperature of 90° C. Dispersions containing 3.5, 6.0, 7.0, 8.0, 
and 11.0% solids gave average loaf volumes of 815, 775, 810, 770, and 775 ce. 
respectively. The results show no appreciable differences between the various 
preparations. It seems probable therefore that solutions containing at least 
11% solids can be dried satisfactorily. Although concentrations greater than 
11% were not tested it appears that such dispersions could be dried without 
damage if the atomizing system of the particular drier used was capable of 
producing the necessary uniform fine spray from the somewhat more viscous 
material. 


Baking Characteristics of Spray-dried Gluten 


Experiments were done to evaluate more fully the baking strength of spray- 
dried preparations obtained under optimum drying conditions. Gluten dis- 
persions prepared in the usual way from a straight run flour,* containing 
14.5% protein and 0.042% ash, were dried with an inlet air temperature of 
160° C. and outlet temperature of 90° C. Using a cake flour containing 7.3% 
protein as a base, samples of the A fraction were added to give a series of 
flours ranging from 7.3% to 22% protein (14% moisture basis). Results of 
baking tests on these are given in Fig. 2. The straight line was calculated from 
experimental points by the method of least squares. Application of conventional 
statistical methods gives a regression coefficient, for loaf volume on per cent 
protein, of 48.2 with a standard error of 1.84. The standard deviation of the 
loaf volumes from the regression is 32.2 cc. 

The regression coefficient of 48.2 indicates good baking quality and com- 
pares favorably with results of Aitken and Geddes (1), who, using dried glutens, 
reported values ranging from about 30 to 59 depending upon the baking 
formula used. Although Aitken and Geddes obtained an average coefficient of 
30.7 without addition of bromate, preliminary experiments in the present 
investigation indicated that no improvement was obtained upon addition of 
oxidizing or reducing agents. 


*Supplied by Robin Hood Mills, Saskatoon, Sask. 
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Fic. 2. The effect on loaf volumes of adding spray-dried gluten to a 7.3% protein cake 
flour. (Regression coefficient for loaf volume on per cent protein = 48.2+1.8. Standard error 
of loaf volume 32 cc.) 


Comparison of Spray-dried Gluten with a Commercial Vital Gluten 

The baking behavior of spray-dried gluten was compared with that of a 
commercially available gluten (hereinafter called vital gluten) of the type 
used in some dietetic breads. Results with cake flour, raised arbitrarily to 15% 
and to 40% protein levels by supplementation with the gluten preparations, 
are given in Table II. 

At the 15% protein level results with spray-dried gluten are similar to 
those presented earlier; however loaf volumes with “‘vital gluten” were below 
600 cc. The loaves were hard and stiff textured and would be quite unacceptable 
for human consumption. At the 40% protein level, ‘‘vital gluten” gave the 
most acceptable loaves. Those from spray-dried gluten were open textured, 
soft and gummy, and had volumes greater than 2000 cc. 


TABLE II 
COMPARISON OF FLOURS SUPPLEMENTED AT TWO LEVELS WITH DRIED GLUTENS 





Protein Average 
content, % loaf 
Composition of flours (14% volume, Comments 
moisture) cc. 
Cake flour only 7.3 420 
Cake flour + “vital gluten” 15.0 580 Stiff texture, gray color 
Cake flour + spray-dried gluten 15.0 800 Even texture, good loaf 
Cake flour + ‘‘vital gluten”’ 40.0 710 Smooth texture, brownish color 
Cake flour + spray-dried gluten 40.0 Over 2000 Very open texture, soft, and 


gummy under pressure 
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Although spray-dried gluten evidently is not suitable for preparation of 
protein-rich bread, the results shown in Table II are another demonstration 
of its baking strength. Using the volumes at 7.3% protein and 15% protein, a 
regression coefficient of 49 is estimated. (See Fig. 2.) If the linear relation 
between protein content and loaf volume holds over a wide range of concen- 
trations, the volume expected at 40% protein would be 2020 cc. This is con- 
sistent with the estimated volumes above 2000 (Table I1). 


Water Vapor Adsorption by Dried Gluten } 


In connection with problems of packaging and storage, water vapor ad- 
sorption by gluten at different humidities was investigated using apparatus as 
described by Colvin, Craig, and Sallans (2). Humidifying columns thermo- | 
stated at 25° C. were used to condition the atmosphere for 1-gm. samples of 
previously dried gluten. Air of controlled humidity was passed over the samples | 
at the rate of 360 ml. per hour until moisture equilibrium was obtained (indi- 
cated by constant weight). Four to five days were usually required for this. | 

The results in Table III show a pattern of moisture adsorption typical of 
plant materials (2) and suggest that dried gluten will adsorb some moisture 








TABLE III 
EQUILIBRIUM MOISTURE CONTENT OF GLUTEN AT VARIOUS 
RELATIVE HUMIDITIES (% BY WT.) | 


Equilibrium moisture content, % 






























Relative Saturated salt —_—_—_—___———_ — —_—————— 

humidity solution in Spray-dried gluten Freeze-dried gluten 
of air, humidifying — = ~- + - 
% column Sample 1 Sample 2 Sample 1 Sample 2 

11 LiCl 5.5 4.0 4.5 3.9 

33 MgCl. 6.9 7.4 5.4 7.6 

52 Mg(NOs)2 12.0 10.0 9.4 1 

75 NaCl 13.4 14.7 12.9 17.1 

93 KNO; 17.9 18.8 18.7 16.5 








from the atmosphere even when stored at relative humidities as low as about 
10%. At humidities of 93%, the samples became wet in appearance and tended 
to shrink to a gummy mass. The uptake of moisture even at fairly high hu- 
midities is however not excessive and under ordinary conditions would not 
likely give rise to serious difficulties in storage and packaging. 








Microscopic Examination of Spray-dried Gluten 





To further characterize the spray-dried glutens, photomicrographs of a 
number of preparations were made. From these, estimates of the particle size 
and size distribution were obtained. Fig. 3 shows typical photomicrographs of 
spray-dried gluten and of household flour. Photograph A shows that the 
powder particles are spherical in shape and not greatly variant in size. A 
number of particles appeared as hollow spheres with ‘‘vapor”’ holes in their 
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Fic. 3. Photomicrographs of some spray-dried glutens and of household flour. 
A—spray-dried gluten; B—household flour. 


side (see arrow). The hollow spheres are characteristic of many spray-dried 
materials (5). 

An estimate of the particle size distribution in some spray-dried gluten 
preparations was obtained by scaling 100 randomly selected diameters on 
photomicrographs. Between 50 and 60% of the particles were between 10 and 
20 microns in diameter. In most samples 20% or more of the particles were 
less than 10 microns in diameter and although there were always some of 
larger size, there were very few with diameter greater than 40 microns. Photo- 
micrograph A was chosen to show the appearance of the particles rather than 
their size distribution. 

The relative fineness of the gluten suggests that it should be uniformly 
distributed throughout the flour by dry mixing (compare photomicrograph B), 
and indicates that little would be gained, as far as uniform reincorporation of 
gluten is concerned, by smaller particle size. 
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Summary of Properties of Spray-dried Gluten 

Table IV is a summary of characteristics of spray-dried gluten. The almost 
complete dispersibility in acetic acid is in accord with evidence that the 
physical characteristics of gluten are not extensively altered by its separation 
and spray drying. 


TABLE IV 


SOME PROPERTIES OF SPRAY-DRIED GLUTEN 


Nitrogen content 15.4%+0.3% 
Moisture content as prepared 0.4% to 4% 
Equilibrium moisture content: 

Relative humidity 33% 7% 

Relative humidity 75% 14% 
Bulk density (flour 0.6) 0.17+0.03 


Per cent of total nitrogen dispersible in 0.01 N acetic acid (8% solids) 99.4 
Regression coefficient: 

Loaf volume over % gluten added 48 
Median particle diameter (av.) 18 microns 





In summary, results presented herein provide evidence that wheat gluten 
retaining good baking quality can be prepared by spray drying acetic acid 
dispersions of the protein. Since a good product was obtained over a range of 
drying conditions, it would appear that, with the more efficient large-scale 
drying units, satisfactory control of drying variables could be readily achieved. 
The experimental data presented on several physical characteristics of spray- 
dried gluten may be of interest when considering the utilization of the material 
or the problems of its packaging and storage. 
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HUMIDIFICATION OF COLD STORAGES: THE JACKET SYSTEM! 


By C. P. LENTz 


ABSTRACT 


Basic principles of the jacket system of maintaining high relative humidity in 
cold storages are outlined. Drying by the cooling coil is prevented by circulation 
of the refrigerating air stream through a jacket surrounding the room rather 
than through the room. Relative humidities between 89 and 98% were main- 
tained in the test room, the level depending mainly on rate and direction of 
air circulation in the jacket and on air flow resistance. Spatial variations in room 
temperature were less than 2° F. during most tests and velocity of air movement 
was negligible. Reduced frosting of the cooling coils in jacketed systems may 
offset most of the 10-15% increase in operating cost resulting from the increased 
air circulation required. The application of the test results in the design of 
jacketed storage rooms is discussed. 


INTRODUCTION 


The use of a jacketed storage room was first proposed by A. G. Huntsman 
in 1931 as a means of reducing surface desiccation of frozen fish (2). Tests 
reported by Huntsman showed that drying in the room he designed was 
reduced to one-fifth of that occurring in conventional storage rooms. Further 
work on the jacketed system was reported by O. C. Young in 1933 (5). His 
work showed that uniform temperatures could be maintained in a jacketed 
room, and emphasized the importance of an effective vapor seal between the 
room and the jacket air space. This type of room, with its uniform temperature 
and high relative humidity, would appear to have a marked advantage over 
conventional rooms in the storage of most frozen foods. Although a number 
of jacketed storages have been built, the principle has not been applied ex- 
tensively, presumably because of the lack of generally applicable design and 
operating information. This paper reports an investigation made to obtain 
such information. 


Principle of Operation of the Jacket System 


Surface desiccation is caused mainly by the migration of moisture from the 
surface of the stored product to the refrigerating surface, the rate of desiccation 
being proportional to the difference in vapor pressure between produce surface 
and refrigerating surface. In the jacketed system the refrigerating air circu- 
lates through a jacket or system of ducts entirely surrounding the room, 
rather than through the room itself (Fig. 1). Hence, the difference in vapor 
pressure between product and refrigerating surfaces is not directly dependent 
on the external heat load or the operation of the refrigeration system. 

If no heat or moisture is given off or absorbed in the room and if temperatures 
in the jacket remain constant, the equation 


Q = WC(tze — hh), 


1Manuscript received March 24, 1958. 
Contribution from the Division of Applied Biology, National Research Laboratories, Ottawa. 
Issued as N.R.C. No. 3630. 
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Fic. 1. Schematic sketch of jacketed room. 


where Q = heat transferred into room through insulation, 
W = weight of air circulated through jacket, 
C = specific heat of air, 
t; = temperature of air entering jacket, 

tp = temperature of air leaving jacket, 


can be used to explain the operation of the system. If air flow distribution in 
the jacket is uniform, the room temperature will be (¢; + t2)/2 and if the 
inside wall of the jacket is vapor-proof but without thermal resistance the 
lowest surface temperature and dewpoint of the air in the room will be h. 
For a given Q, it is apparent that the relative humidity in the room will depend 
on W, the amount of air circulated through the jacket. 

In practice, the relative humidity in the room may be affected by variations 
in the temperature of the air in the jacket, uneven distribution of the circu- 
lating air, and heat loads in the room. Since these factors determine the amount 
of air that it is practical to circulate in the jacket, they will be considered 
separately. 


(a) Temperature Variations in Jacket Air 

Temperature variations in the jacket tend to reduce the humidity in the 
room, The magnitude of the effect depends on the amplitude and period of the 
variations, and on the thermal inertia of the inner wall of the jacket. As 
pointed out by Young (5) the use of some thermal insulation in the jacket wall 
(e.g. 2 in. wood sheathing) will reduce the effect of temperature cycling of the 
jacket air. 


(b) Heat Loads in the Room 

The effect of heat loads in the room is to raise the room temperature above 
(t; + t)/2, and hence to lower the relative humidity. Heat transfer into the 
room from the jacket has a similar effect. Since there is a temperature gradient 
in the jacket air (t to fa), heat will flow into the room through the warmer 
surfaces, and out through the colder surfaces. If the inner jacket wall has no 
thermal resistance, the coldest part of the wall surface will be ¢, and the 
warmest part fa. Insulation will reduce the heat flow and hence lower the 
warmest surface temperature and raise the coldest surface temperature. 
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Heat transfer into the room through the structural members in the jacket air 
space is also reduced by this insulation. 


(c) Air Flow Distribution in Jacket 

Uneven air flow distribution in the jacket will act in two ways to reduce the 
relative humidity. First, additional heat will be transferred into the room from 
parts of the jacket having reduced flow, and second, the area through which 
heat is transferred out of the room and into the jacket will be reduced. Both 
of these factors tend to increase the difference between the room temperature 
and the temperature of the coldest part of the jacket wall, and hence to lower 
the relative humidity. It can be seen that if temperature variations and uneven 
flow distribution in the jacket result in a room temperature 5° F. above the 
coldest wall surface temperature, there is slight advantage in circulating 
enough air to reduce (f2 — t;) to 1° F. 


APPARATUS AND PROCEDURE 
Description of Test Room 
The test room (Fig. 2) was constructed so that the results obtained could 





Fic. 2. Cut-away view of jacketed test room. 
(1) Vapor barrier. 
(2) Mineral wool insulation 8 in. thick. 
(3) Ceiling air space. 
(4) Wall ducts. 
(5) Floor air space. 
(6) Insulated and vapor proof inside wall (asphalt painted plywood over ‘‘Ten-test 
3 in. thick). 
(7) Ceiling and floor air guides, 6 in. deep. 
(8) Wall air guides, 23 in. deep. 
(9) Top plenum chamber. 
(10) Bottom plenum chamber. 
(11) Horizontal partition. 
(12) Dampers to change direction of air circulation. 
(13) Air in. 
(14) Air out. 


” 
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be applied readily in the design of commercial storages. The inside dimensions 
of the room were 10 X 10 X 20 ft. long. The system was designed so that air 
could be circulated through the jacket from top to bottom or from bottom to 
top at any rate desired. The jacket ducts were constructed so that uniform 
flow would be obtained in the walls rather than in the entire jacket. This was 
done by reducing the depth of the ducts in the walls compared with those in 
the floor and ceiling, and thus increasing the proportion of the total friction 
pressure loss occurring in the walls. Since ideal distribution requires the 
quantity of air flowing in each circulation path to be proportional to the 
length of the path, the design used represents a compromise in the interest of 
practical simplicity. The inside wall of the jacket was insulated and vapor- 
proofed for reasons already mentioned. 

Accurate temperature control was obtained by maintaining a constant 
ambient temperature and by heating the air entering the jacket to eliminate 
the effects of changes in refrigeration capacity. 


Measurements 

Temperatures were measured with thermocouples installed so that con- 
duction errors would be minimized. Twenty-seven elements were located in 
the room air, seven in the jacket air, and ten were attached to the inside surface 
of the room. The surface thermocouples were placed in areas of maximum 
and minimum surface temperature. In addition, fourteen differentially con- 
nected couples were used to measure temperature differences between top 
and bottom of wall ducts. Variations in flow distribution in the ducts were 
estimated from these differences. 

Relative humidity was measured at five locations in the room by the wet- 
and dry-bulb method. This method had been found to give satisfactory results 
in an earlier investigation (3). Each instrument consisted of a dry-bulb 
thermocouple and a series of ten differentially connected wet- and dry-bulb 
junctions. The measurements were checked by the gravimetric method using 
chemical absorption. Relative humidity was calculated with reference to the 
vapor pressure of ice since this is the datum of interest in frozen food storage. 

Velocity and static pressure measurements were made with a micromano- 
meter and Pitot tube. Air flow rates were calculated from velocity pressure 
measurements made in an orifice in the inlet to the jacket plenum chamber 
(Fig. 2). Jacket friction pressure loss was calculated from the pressure differ- 
ence between the top and bottom sections of the plenum chamber. 


Scope of Tests 

The effects of jacket air flow rate, duct friction pressure loss, and forced air 
circulation in the room were studied with both top-to-bottom and bottom-to- 
top circulation in the jacket. The rate of air flow in the jacket was varied to 
give differences in temperatures of 1 to 10° F. between the air entering and 
leaving the jacket. This range includes the entire range likely to be of interest 
in storage practice. The effect of duct resistance was studied after preliminary 
tests indicated the importance of this variable. This resistance was varied by 
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partially closing the dampers in the wall ducts. Tests were also made using 
forced circulation in the room because such circulation may be necessary in 
practice when produce must be cooled quickly. The amount of circulation 
during these tests was approximately equal to that in comparable rooms using 
forced air cooling units. In addition, two tests‘were made to determine the 
value of insulating the jacket wall. In one of these, the temperature of the 
air entering the jacket was varied through a +4° F. cycle in four hours and in 
the other a heat source equal to half the external heat leakage was put in 
the room. 

The room was empty and at a temperature of 0° F. for all of the tests. The 
room was not loaded because the thermal inertia of a load would have de- 
creased the sensitivity of the room to changes in operating conditions. Other 
temperatures were not tested because the performance of the room at one 
temperature can be used to predict its performance at other temperatures. 

Most of the tests were continued for two or three days after constant 
temperatures were reached. Several longer tests (two weeks) were made to 
ascertain whether or not any slow changes were taking place. In all, the 
jacketed room was operated for more than a year. 


RESULTS 
Relative Humidity 

Relative humidity in the room varied with rate and direction of air circu- 
lation in the jacket, decreasing from 98% at high rates of circulation to 89% 
with a low rate of flow in the top-to-bottom direction (Fig. 3). The decrease 
was gradual with bottom-to-top circulation and did not exceed 2% at the 
lowest rate of flow tested. With top-to-bottom circulation the relative hu- 
midity dropped sharply as the flow rate was reduced below 1200 c.f.m. This 
critical circulation rate was reduced by partial closure of the wall ducts with 
dampers. The dampers had little effect with bottom-to-top circulation. 

Since the dewpoint in the room was higher than the lowest wall surface 
temperature during most of the tests, relative humidities based on minimum 
surface temperature (Fig. 4) are lower than the actual relative humidities 
measured (Fig. 3). This discrepancy was caused by the slow migration (from 
warmer to colder wall surfaces) of moisture deposited on all surfaces during 
door openings. Since the dewpoint temperature in the room would eventually 
equal the lowest surface temperature, the calculated relative humidities are a 
better criterion of performance than those obtained by the wet-bulb method. 
Top-to-bottom circulation gave the best results on the basis of minimum 
surface temperature. 


Air Movement in the Room 

Smoke tests showed that air velocities in the jacketed room were negligible 
(less than 5 f.p.m.) under all conditions tested. 

Forced air circulation in the room reduced relative humidities calculated 
from wet-bulb measurements (Fig. 3) but had no appreciable effect on those 
based on minimum surface temperature (Fig. 4). Since earlier experiments (3) 
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Note: (1) For approximate air velocity in wall ducts divide c.f.m. by 10 (Figs. 3-7 incl.). 

(2) ‘Temperature difference’”’ refers to increase in temperature of air in circulating 
through jacket (Figs. 3-7 incl.). 


showed that similar changes in air velocity (5 f.p.m. to 100-300 f.p.m.) did 
not cause error in the wet-bulb depression, the results indicate that forced air 
circulation decreased the dewpoint of the room air. The decrease may have 
been caused by a proportionally greater increase in the rate of condensation 
than in the rate of evaporation of migrating moisture. 


Spatial Temperature Variations 

Spatial variations in room temperature were larger with bottom-to-top 
circulation than with top-to-bottom circulation (Fig. 5). The critical rate of air 
flow indicated by the relative humidity measurements is also shown by the 
temperature variations. Variations in dewpoint temperature (calculated from 
wet-bulb measurements) were similar to those in the dry-bulb temperatures. 
The dewpoint variations are of no consequence in assessing the performance 
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of the room, however, since the door openings causing these variations would 
not have appreciable effect in practice. 


Jacket Temperatures and Friction Pressure Loss 


The average jacket temperature was approximately the same for top-to- 
bottom and bottom-to-top circulation at high rates of flow (Fig. 6). As the 
rate of flow with top-to-bottom circulation was decreased below the apparently 
critical rate, the average jacket temperature dropped sharply. Temperature 
differences in the wall ducts indicated that the change in operating character- 
istics was due to non-uniform air flow distribution in the ducts. At the critical 
flow rate the direction of circulation in some of the ducts apparently reversed 
and the number of ducts with reversed circulation increased as the rate of 
circulation was further reduced. Calculations based on duct air temperatures 
showed that differences in density due to differences in temperature between 
ducts were sufficient to reverse the direction of circulation at these low rates 
of flow. Increases of 10-15° F. in surface temperatures in some locations gave 
further evidence of poor circulation. It may be noted that the difference 











CANADIAN JOURNAL OF TECHNOLOGY. VOL. 33 


- 2.0 

Ww 

a 

=» 

-a 

< ~TOP-T0- BOTTOM 

xe 10 { 

Ww 

a a at ist li Ni 

: Pee 

Ww 

- oO ss ’ 

w 

° og F777 0--0- OO 
a i 

= +10) o Z gottom-to-toP 

> 

<a 


BOTTOM 


(8) 


‘gS ee BOTTOM-TO-TOP 
2.0 Qe. 


“ © 
Oe wee --1-0.-- 


1.0 


° 


\—TOP-TO-BOTTOM 


400 800 1200 1600 2000 2400 
JACKET AIR FLOW RATE, C.F.M. 


SPATIAL TEMPERATURE VARIATION,°F 


15 10 5.0 40 3.0 2.0 1.5 
TEMPERATURE DIFFERENCE, °F 


Fic. 5. Spatial temperature variation in test room. (A) vertical, (B) maximum. 
@ © Ducts open. §§ 0 Ducts partially closed. . 
Closed symbols: top-to-bottom circulation. Open symbols: bottom-to-top circulation. 


between the average jacket temperature (¢; + ¢)/2 and the average room 
temperature (0° F.) is a measure of the amount of heat flowing into the room 
through the jacket and its structural members. 

Friction pressure loss in the jacket was doubled (Fig. 7) by partial closure 
of the wall duct dampers to improve air flow distribution in the jacket. Even 
at high rates of flow, however, the friction pressure loss in the jacket was not 
significant when compared with the pressure loss in the evaporator and fan 
(0.5 in water). 


Temperature Cycling and Heat Load Tests 

During the temperature cycling tests the room temperature varied one-third 
as much as the jacket temperature. It should be noted that the variation in 
room temperature would be even smaller for cycles of shorter period. 

A heat load of 2000 b.t.u./hr. raised the room temperature 5.7° F. This is 
approximately the amount the same heat load would raise the temperature 
of a comparable room cooled by a forced air circulation evaporator operating 
continuously with a 10° F. refrigerant-room air temperature difference. 
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DISCUSSION 

Relation of Operating Characteristics to Jacket Air Circulation 

Relative humidity, minimum surface temperature, and spatial temperature 
variations in the room depended on the rate and direction of air circulation 
in the jacket. With bottom-to-top circulation these variables changed at an 
approximately uniform rate as the rate of air flow was changed. With top-to- 
bottom circulation there appeared to be a critical rate of flow below which 
changes occurred more rapidly. The critical rate of flow appeared to depend 
on flow resistance. These characteristics of the room can be explained by 
considering the effect of two factors on the operation of the system. These 
factors are the constant rate of heat flow into the jacket air space, and the 
differences in flow resistance between air flow paths in the jacket. Less air 
flows in those paths having higher resistance than others. Since the same 
amount of heat is absorbed in each of the wall ducts, the average temperature 
in the ducts with lower rates of flow will be higher than in others. With bottom- 
to-top circulation, the resulting differences in density between ducts tend to 
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maintain uniform flow. With top-to-bottom circulation, however, these dif- 
ferences in density tend to increase differences in flow. When the pressures due 
to differences in density approach the friction pressure loss (Fig. 7), flow 
becomes unsteady and may reverse in direction. The resulting poor distri- 
bution of refrigerating air in the jacket increases spatial temperature variations 
in the jacket and the room, and hence, lowers relative humidity. 


Design of Jacketed Rooms 

If allowance is made for changes in dimensions, insulation, and temperature 
conditions, the results obtained with the test room should be generally applica- 
ble in the design of jacketed rooms. Although the jacketed room represented 
an improvement over most directly cooled rooms with all of the types of 
operation tested, the results showed that direction and rate of air circulation, 
duct design, and jacket wall construction affected conditions in the room 
significantly. Each of these factors will be considered separately. 


(a) Direction and Rate of Air Circulation 

With good distribution of flow in the ducts, top-to-bottom circulation gave 
higher minimum surface temperatures and smaller spatial temperature varia- 
tions than bottom-to-top circulation (Figs. 4 and 5). The resulting difference 
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in rate of surface desiccation would be increased by forced circulation in the 
room. It may be noted that the advantage of top-to-bottom circulation de- 
pends on uniform circulation in the jacket. With poor construction, or insuffi- 
cient circulation, bottom-to-top circulation could give better results. 

The rate of air circulation required depends mainly on the conditions desired 
and on the heat load to be absorbed. The results indicate that 1000 c.f.m. 
would be the minimum satisfactory rate in the test room and that little would 
be gained by increasing the circulation rate above 1800 c.f.m. These limits 
correspond to temperature differences of 4 and 2° F. respectively, between the 
air entering and leaving the jacket. The cost of fan and refrigeration power 
should also be considered. Since lower jacket temperatures require lower 
refrigerant temperatures, the cost of refrigeration power increases as the rate 
of circulation is reduced. Fig. 8 shows that the cost of fan power increases 
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Fic. 8. Relation between jacket temperature difference and air circulation rate. 


rapidly as the jacket temperature difference is reduced below 2° F., and the 
refrigeration power increases rapidly as the rate of air circulation is reduced 
below 1200 c.f.m. per ton, on a sensible heat basis. 


(b) Duct Design 

The main object of the jacket ducts is to distribute air through the jacket 
in proportion to the heat to be absorbed. For top-to-bottom circulation a 
practical design appears to be one in which the friction pressure loss in the 
wall ducts is uniform, represents a large proportion of the total friction 
pressure loss in the jacket, and is appreciably larger than the pressure differ- 
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ences due to differences in density between ducts. Increasing the duct re- 
sistance in the test room improved its operation, and the results indicate that 
further increase in resistance would have given further improvement. The 
maximum resistance that may be used depends only on cost of fan power. 
In the test room the resistance of the jacket at the critical rate of flow was 
less than 0.01 in water. Since this represents less than two per cent of the 
pressure loss normally occurring in evaporator-and-fan systems, it could be 
increased by a factor of 5 to 10 without causing a prohibitive increase in fan 
power. Increased resistance should be obtained by using smaller ducts rather 
than greater flow, because greater flow also increases pressure loss across the 
evaporator. Fig. 9 shows that the test room should have been designed with 
ducts 1-1} in. deep instead of 23 in. deep. 

The minimum permissible flow resistance is less easily decided because it 
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Note: (a) For approximate air velocity in wall ducts 2} in. deep divide c.f.m. by 10, veloc- 
ities in other wall ducts inversely proportional to duct depth. 
(6) Assumptions made in calculations: 
. Jacket and wall duct entrance loss = 0.5 v.p. 
. Jacket and wall duct exit loss = 1.0 v.p. 
. Area of jacket entrance and exit = 5 sq. ft. 
. Cross-sectional area of wall ducts 2} in. deep = 11 sq. ft. 
. Friction pressure loss in floor and ceiling jacket negligible, except for entrance 
and exit losses. 
6. Absolute roughness of jacket surfaces 0.05 in. 
. The basic flow equation for circular ducts is applicable to the circular equiva- 
lents of the wall ducts. 
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depends mainly on the variation of flow resistance in the ducts, which is a 
function of materials and construction, as well as design. Wall height and the 
ratio of air circulated to heat absorbed also affect the minimum pressure loss 
which will give satisfactory operation. A higher rate of heat absorption gives 
larger temperature differences and hence larger differences in density and 
pressure. Although both friction pressure loss and pressure difference due to 
difference in density are proportional to duct height, allowance should be 
made for duct height because of the proportion of the duct resistance repre- 
sented by entrance and exit pressure losses. The minimum permissible loss 
should cause no difficulty in design however, since a liberal factor of safety 
can be allowed without appreciable increase in cost. In addition, friction 
pressure loss can be adjusted with dampers if necessary. 


(c) Jacket Wall Construction 

The main requirements of the jacket wall are that it be a good water vapor 
barrier, and that it have some thermal resistance, as already discussed. The 
amount of insulation which should be used depends on the conditions desired. 
For most purposes the insulation used in the test room would be adequate. 


CONCLUSION 


The results of this study indicate that surface desiccation of frozen products 
during storage can be reduced to negligible levels by use of the jacketed room. 
Not only can relative humidities approaching 100% be maintained, but the 
uniform temperatures and low air velocities in the room should reduce surface 
desiccation to levels much below those which would occur in directly cooled 
rooms at similar levels of relative humidity. 

In addition to reducing surface desiccation of stored produce, the jacketed 
room offers a solution for the problem of deterioration of insulation resulting 
from moisture condensation. This problem is particularly serious where the 
vapor pressure and concentration gradients periodically reverse in direction. 
Moisture condensation in insulation was discussed by Munters (4) who sug- 
gested the circulation of dry air through the insulation. 

The operating cost for a jacketed room may not be much higher than for a 
comparable directly cooled room. Although the increased air circulation re- 
quired for jacketed operation increases power consumption by 10-15%, this 
increase in cost may be more than offset by the increased efficieacy of operation 
of the evaporators that results from reduced frosting. 

In the jacketed system high relative humidity is maintained by preventing 
moisture from leaving the room. Since this method is not effective where there 
is moisture absorption in the room, moisture absorption by packaging ma- 
terials may cause appreciable surface desiccation. Such absorption would not 
be a problem in a system in which high relative humidity is maintained by 
moisture addition (as in the direct humidification system described earlier (1, 
3). The answer to the problem may be provided by a combination of the two 
systems now being studied. 
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AUTOMATIC TRANSLATION OF AERATION DATA BETWEEN 
DIFFERENT FERMENTORS! 


By Pinc SHu 


ABSTRACT 


The present paper describes an apparatus with which the rate of oxygen up- 
take of a fermentation conducted in a master fermentor such as shake flask may 
be automatically duplicated in a slave fermentor, such as a deep tank. The 
device consists of an element coupling the fermentors and a closed loop aeration 
control system for the deep tank fermentor. The function of the former is to 
compare the oxygen uptake rate of the fermentor with that of the shake flask. 
The function of the latter is to correct the difference in the oxygen uptake rates 
transmitted from the comparator by regulating the agitator speed of the tank. 
The agitator speed and the power requirement for such duplication are de- 
termined. Hydrogen peroxide is used as a source of oxygen supply for the fer- 
mentor and gaseous oxygen for the shake flask; the scale of fermentation may be 
amplified by varying the hydrogen peroxide concentration. The performance 
of this translation device was demonstrated with fermentations by Ustilago seae 
and Aspergillus niger cultures. The device offers a solution to the problem of 
translating the aerations in the event of a change in the scale of a fermentation 
or the type of the fermentors employed. It also provides a convenient means for 
direct comparison of the aeration efficiency of different types of fermentors of 
various sizes under actual fermentation conditions. 


INTRODUCTION 


Among the nutrients utilized during an aerobic fermentation only oxygen 
is supplied continuously. The method for controlling its rate of supply varies 
with the type of the fermentor employed. 

In shake flask fermentations, the extent of aeration may be altered by 
changing either the speed of the shaker or the surface volume ratio of the 
medium. In the deep tank fermentor, aeration is varied by changing either 
the stirrer speed or the air flow (2, 9). The relation between these variables and 
the rate of oxygen transfer through the gas-liquid interface is complicated. 
It varies with different fermentations and involves such parameters as surface 
tension, and dissolved and suspended materials. 

A convenient sulphite oxidation method was introduced by Cooper et al. (3) 
to establish the relation between stirrer power and rate of oxygen transfer in an 
aerated tank. This technique has been used by Karow et al. (6) to scale up a 
streptomycin fermentation from 5 liters to 15,000 gal. Wegrich and Shurter (13) 
scaled a penicillin fermentation from 5 gal. to 24,000 gal. by using stirrer 
power data alone. Since these procedures do not involve the measurement of 
the rate of oxygen transfer during actual fermentation, they can only be 
applied to geometrically similar fermentors. Hixson and Gaden (4), Wise (14), 
and Bartholomew et al. (1) employed a polarographic method to determine the 
oxygen transfer rate in presence of actively metabolizing organisms. This 

1Manuscript received March 31, 1956. 

Contribution from the National Research Council of Canada, Prairie Regional Laboratory, 
Saskatoon, Saskatchewan. Issued as N.R.C. No. 8629. Presented at the Biochemical Engi- 


neering Symposium of the 47th Annual Meeting of The American Institute of Chemical Engineers, 
New York, N.Y., December 12-15, 1954. 
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permits the comparison of aeration efficiency of different fermentors under 
actual fermentation conditions. Recently Hoover et al. (5) described a device 
with which the oxygen concentration in the effluent gas from a fermentor 
may be continuously recorded by measuring its paramagnetic susceptibility. 
This offers a more convenient means for studying the aeration problem. 

In the development of a fermentation process, it is often necessary to change 
the size and type of fermentor. For example, shake flasks are commonly used 
in the initial stage of development and a small tank fermentor at the next 
stage. It is a difficult problem to find the conditions that will reproduce the 
magnitude of aeration. 

The present paper introduces a device that permits the rate of oxygen up- 
take of a fermentation in a given fermentor to be automatically duplicated in 
another fermentor of different size and type. The stirring speed required and 
the power demand to obtain such aeration are determined. It provides a 
convenient means of finding the conditions for reproducing the magnitude of 
aeration in the event of enlargement of the scale of a fermentation or a change 
in the type of the fermentor employed. Specific examples are given to illustrate 
the duplication of the oxygen uptake rate of a shake flask fermentation in a 
deep tank fermentor. 


APPARATUS 
General Description 


To translate the aeration data from one fermentor to another, the fermen- 
tation must be done in these two fermentors simultaneously. The rates of 
oxygen uptake are compared from time to time. Whenever differences occur, 
the oxygen uptake in one fermentor is adjusted so that it follows that of the 
other. Conditions to obtain the required magnitude of aeration may then be 
determined. Fig. 1 shows a schematic diagram of such a system in which a 
shake flask se-ves as the master fermentor while the rate of the oxygen uptake 
is duplicated in a small stirred tank. The closed shake flask system (block A) 
is the same as that described in a previous paper (10) except that the position- 
ing motor only moves the level follower of the oxygen reservoir upward. The 
apparatus for measuring the oxygen uptake in the tank fermentor (blocks B 
and D) is similar except that oxygen is supplied by the catalytic decomposition 
of hydrogen peroxide. Lead, lead oxide, permanganate, and selenium in fused 
granules are used as catalysts (7). The decomposition of hydrogen peroxide is 
instantaneous and quantitative. The amount of water collected in the cylinder 
is a measure of the volume of hydrogen peroxide consumed. 

Air in the fermentor is recycled with a diaphragm pump? after it is passed 
through 30% sodium hydrozide solution. Part of the air is recycled over the 
top of the fermentor and the rest passes through the medium. The extent of 
total recycling is adjusted with valve V2 on the pass shunt across the pump. 
The rate of air passing through the medium is varied by adjustment of valves 
V, and V3. 


2Dyna-pump, a product of Scientific Glass Apparatus Co., Bloomfield, N.J. 
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During a fermentation, the oxygen consumed by the organisms upsets the 
pressure equilibrium between the fermentor and the gas spaces of water 
collector and oxygen generator. An equivalent amount of hydrogen peroxide 
is hence drawn in to compensate for the quantity of oxygen utilized. The 
concentration of hydrogen peroxide is selected to compensate for the difference 
in volume between fermentor and shake flask. 

As the water levels in both the oxygen reservoir and the water collector 
cylinders rise, their rates of rising are compared from time to time with the 
help of a level comparator. The difference is amplified mechanically to operate 
relay 3, which with the help of a positioning motor 3* turns the variac speed 
control dial of motor A in such a direction that it changes the stirring speed to 
nullify the difference. A synchronizing timer switch is used to close the circuit 
between relay 3 and positioning motor 3 for a preset length of time after either 
one of the level follower bulbs makes a move. The speed of motor A is coupled 
to that of the stirring motor B through a speed comparator and a relay which 
operates the positioning motor 4 for the speed control dial of the stirring motor 
(block C). This arrangement eliminates the effect of the load disturbances on 
the stirring speed. The latter responds only to the command of the set speed 
of the motor A. Thus, the speed of stirring in the fermentor is indirectly 
regulated by the relative rates of the level rise in the cylinders. Changes in the 
stirring speed alter the oxygen uptake rate of the organism in the fermentor. 
This rate is then expressed in terms of the rate of rise of the level in the water 
collector cylinder. This closed loop regulating system makes it possible for the 
consumption of hydrogen peroxide by the tank fermentor to follow the oxygen 
uptake in the shake flask. 


Level Comparator and Error Amplifier 

It may be seen in Fig. 2 that the device to compare water levels in the 
reservoirs and amplify their differences is simply a lever system. The extension 
of the platform which carries the level follower bulb of the water collector 
supplies the fulcrum for the lever arm C. The latter is hinged to the platform 
which carries the follower bulb of the oxygen reservoir. A plastic idler bar B 
which has whisker and main contacts is hinged to the fulcrum rod and presses 
against a friction board with a sponge rubber pad in between. This makes it 
possible for the idler bar to stay where it is left. The lever arm is made of a 
brass bar with a } in. slot cut at the end. The amplification factor is about 3. 
The main contact of the idler bar is placed in the slot. As one of the level 
follower bulbs makes a move, the slot end of the lever bar makes an amplified 
motion. If the amplitude of the motion is beyond the limit of the slot, the lever 
bar will touch the main contact of the idler bar and push the latter to a new 
position. The whisker contact only touches the lever bar when the idler bar 
assumes a lower position. The lever bar is electrically connected to contact B, 
while the main and whisker contacts are connected to contacts A and C re- 
spectively, of the electronic circuit (Fig. 5) which operates the positioning 


‘One revolution per minute Bodine K.Y.C. 22 motor, a product of Bodine Electric Co., Chicago, 
Til. 
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motor 3. The connections are made such that when the idler bar assumes 
lower position the positioning motor turns the variac dial to increase the 


stirring speed and vice versa. 
The variac dial is turned by the positioning motor 3 with friction drive. 
The ratio of the diameter of the friction disks is 1:17. The diameter of the 
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variac dial disk is 8.5 in. A 36 gage nichrome wire wound around the edge of 
the dial disk and a stationary sliding contact serve as a potentiometer to record 
the position of the dial. A current passing through the potentiometer is so 
adjusted that the potential drop across the entire wire is 9.7 mv. Since motor A 
carries only a small constant load, the motor speed can be expressed in terms 
of the dial position. 


Synchronized Contact Time Limiter 


When the idler bar assumes a new position the positioning motor will turn 
until one of the level followers moves to disengage the relay contact. Since the 
time interval between two successive movements of the level follower bulbs 
varies widely with the rate of oxygen uptake, this inevitably creates oscillations 
of high amplitude. This variation is minimized by insertion of a switch which 
permits current to flow to the positioning motor for a desired period after 
either or both of the follower bulbs make a move. This sequential process is 
performed by a set of relay arrangements shown in Fig. 3. As mentioned 
previously, the positioning motor relay circuits for the level followers move 
the bulbs only in the upward direction as the water levels rise in the cylinders. 
The other contacts of the single pole double throw relays A and B are available 
to operate normally closed a-c. relays C and D. 

The circuit closure of either or both of these relays (C and D) energizes an- 
other normally closed a-c. relay E, and at the same time opens a solenoid valve 
to fill up the timer cylinder with water, thus closing the circuit at F. The 
relay contact £ and the mercury contact F are connected in series to the power 
line of the positioning motor 3. 

As the positioning motor 1 and/or 2 ceases to move, both the relays C and D 
become energized, the contact E is closed and likewise the solenoid valve. 
The contact F remains closed until water in the timer cylinder drains through 
the overflow 2 to its breaking point. 

The over-all effect of this arrangement is that the positioning motor 3 alters 
the stirring speed only when the main contact of the level comparator, the 
relay contact £, and the mercury contact F make simultaneously closed circuit. 
This can occur only after the level follower bulb makes a move. 


Motor Speed Comparator 


The load on the stirring motor often varies as the physical consistency of 
the fermenting broth changes and causes variations in stirring speed. These 
may be nullified by standardizing the speed of the motor continuously against 
that of another motor running with a small load. This is done by the speed 
comparator as shown in Fig. 4. 

The speed comparator consists of a plexiglass disk A, and a brass disk B 
(each 3 in. diameter). Disk A is driven by a variable speed d-c. motor A‘ after 
a speed reduction with a gear of 1: 800 ratio. Disk B is driven by the d-c. 
stirrer motor B®, with a gear reduction of 1:75. These two disks are mounted 


‘Variable speed motor 2T60, purchased from Gerald K. Heller Co., Baltimore, Maryland. 
51/15 h.p., variac speed control motor, purchased from General Radio Co., Cambridge, Mass. 





wd 
eo 
N 


DATA 


AERATION 


SHU: 


“10jvIeduWIOD paads 10}0JQ“f “DIY 


fmm | 
) 


need fin SY 








@ 9510 Ssvus 


LIVLNOD 
UINSIHA 


42078 Yz70! 


V 9810 Ssv19Ix31d 


sitive 
ONILINGNOD 





” 
L,-) Vora 


SLOVLNOD AVI3Y OL 


v x08 uv39 


‘p puke ¢ s10jOW Bul 


Uoose Av13y 
O1ZW4WNYS Y3LL0d = 2‘ 1 AVT3Y 


vz 3sn4=: 4 
Z2a-'A YO4 SY3LV3H = 2H‘'H 


‘OO AOSI IMB = 2059 
49S = Sy 
Mivoosi = "y'?y 
‘bawi= fy‘'y 
roznzi= 2at'a 








nyisod pue Arjas 





94} 4JOj Wi 


2 AVI3y ——— 


eidseIp INI 


‘OV AGI 





"g ‘DIY 


| AV134y 





286 CANADIAN JOURNAL OF TECHNOLOGY. VOL. 33 


concentrically and turn in the same direction. An idler brass disk (3 in. thick 
and 2.5 in. in diameter) with a contact slot on the side is placed concentrically 
on disk B, and is free to slide. Disk B is connected to terminal B of the relay- 
motor circuit (see Fig. 5). Disk A is furnished with main and whisker electrical 
contacts which are connected to terminals A and C of the relay-motor circuit 
through platinum rails and sliding contacts on the disk. The main contact is 
placed in the slot of the idler block. When disk A rotates faster than B, both 
the main and the whisker contacts are in touch with the idler block and 
activate positioning motor 4 to increase the speed of stirring motor B. When 
the reverse is true, only the main contact touches the idler disk and the posi- 
tioning motor reduces the stirring speed. Thus, the speed of the motor A 
controls that of the stirring motor. 

The coupling between the positioning motor and the variac dial of stirring 
motor is also made with friction disks. To slow down the dial movement, a 
diameter ratio of 1:8.5 was used. The circuit diagram for the relays and 
positioning motors 3 and 4 is shown in Fig. 5. It may be seen that the make 
and break of terminals C and B (whisker contact) reverses the direction of the 
motor while that of the terminals A and B (main contact) switches on and off 
the power line. 


Air Flow Regulator 

During a fermentation, because of changes in stirring speed and increases 
in suspended matter, there is an appreciable variation in resistance to air flow 
through the medium. The resistance variation in air flow is automatically 
corrected with the regulator system shown in Figs. 6 and 7. This regulator is, 
in effect, a capillary flowmeter with a double arm on one side of the capillary 
and a single arm on the other side. The inner walls of these arms are coated 
with water repellent Desicote.* Nichrome electrodes are adjusted to give the 
desired flow rate. They are electrically connected to the contact terminals A, 
B, and C of the relay-positioning motor circuit (Fig. 5). When changes in the 
rate of air flow exceed a previously set limit, the flowmeter fluid in the arms 
will make a contact for electrode B with electrode A in the single arm or with 
electrodes A and C in the double arm. This activates positioning motor 5 to 
turn the air flow control valve V3 in such a direction that it nullifies the effect 
of the disturbance. Air flow is regulated by squeezing the rubber tubing against 
a wall B (Fig. 7). The bolt C (28 threads per inch) is mounted through a 
threaded collar D which is turned in a bearing mounted on a stationary plate 
at a rate of 4 r.p.m. by a 1 r.p.m. positioning motor coupled through a gear 
system G, and Gs. With the rectangular bolt-head A pressing against guide 
plate E, the bolt moves back and forth without turning when the threaded 
collar rotates in one direction or the other. The range of this motion is limited 
by a microswitch pressing against bolt head A. 


Fermentor and Torque Recording Device 
The fermentor, a glass jar (8 liter), is similar to that described by Rivette 
et al. (8) except that the agitator is a four-bladed turbine impeller and the 
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Fic. 6. Air flow regulator. 
Fic. 7. Air flow control valve. 


bearing for the stirring shaft is mounted on a shelf above the fermentor (Fig. 8). 
The entry of the shaft to the fermentor is doubly sealed with mercury and 
lard oil. All tube connections to the fermentor are made with flexible 3/16 in. 
I.D. rubber tubing and are placed close to the center of the jar top. 

The torque measuring device (Fig. 9) is similar to that described by Cooper 
et al. (3) except that the counter weight is a piece of uniform steel rod (9 in. 
long and 0.75 in. in diameter), with added weight on the top, immersed in a 
pool of mercury. The length of rod above the mercury surface represents the 
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magnitude of counter weight required. The change in the extent of immersion 
is measured and recorded as the potential drop across a nichrome wire (28 ohms 
per foot) stretched along the side of the rod. Both wire and mercury must 


be kept clean. 
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EXPERIMENTAL AND RESULTS 


The performance of this device was tested by fermentation of glucose with 
Ustilago zeae PRL 1306 and production of a-amylase with Aspergillus niger 
PRL 558. Both the stability and the reproducibility of the apparatus were 
studied. 


Stability and Reproducibility 


With the present arrangement, oscillation of the stirrer speed may be reduced 
during the period of steady fermentation by decreasing the contact time of the 
synchronized contact time limiter. However, a decrease in contact time would 
hamper the ability of the fermentor to respond to fast changes of the rate of 
oxygen uptake in the shake flask. Furthermore, since the relation between the 
control dial position of motor A and the stirrer speed is non-linear, the effect 
of contact time to the stirrer speed oscillation also varies with the stirrer speed. 
Thus, the choice of a proper contact time is essential. 

With the present control elements, a contact time of five seconds was found 
satisfactory for fermentations requiring stirring speeds higher than 300 r.p.m. 
during the steady fermentation period. 

Data summarized in Fig. 10 shows the stability of stirring speed and torque 
developed during a steady oxygen uptake period of U. zeae fermentation. The 
stirrer speed oscillation was about 10 r.p.m. and the change in torque was 
insignificant. 

The reproducibility of the device was tested with the sulphite oxidation 
method at two levels of aeration. For a given aeration rate, the stirring speed 
may be reproduced within the range of the stable oscillation. 


Glucose Fermentation by Ustilago zeae PRL 1306 


Preliminary runs showed that the initial lag of the fermentation in the shake 
flask was less than that in the fermentor. It was not possible for the uptake of 
oxygen in the fermentor to follow that in the flask even at high rates of stirring. 
Therefore, a culture grown for 24 hr. was used for the coupled fermentation. 

The medium of Thorn and Haskins (12) was used. The flask (500 ml.) 
contained 100 ml. medium and was incubated on a New Brunswick Shaker’ 
(320 r.p.m.) at 26.5° C. The fermentor contained 4 liters of medium, and a 
solution containing 11.8% hydrogen peroxide was employed as the oxygen 
source. Consumption of one volume of this solution by 4 liters of medium is 
equivalent to the utilization of the same volume of gaseous oxygen per 100 ml. 
of medium under these conditions. The initial speed of stirring was set at 
about 100 r.p.m. to ensure that the culture was in suspension. The fermentor 
was kept at 25.5° C. in a water bath. The rates at which air was recycled 
through the medium and sodium hydroxide scrubber were adjusted to 1 liter 
per minute. 

The results of this fermentation are summarized in Fig. 11. When the system 
reached a steady state the stirring speed remained between 325 and 350 r.p.m. 
to duplicate the rate of oxygen uptake in the shake flask (average 200 ml. per 


1A product of New Brunswick Scientific Co., New Brunswick, N.J. 
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liter per hour). During this period, little change in torque occurred and the 
variation of stirring power was caused mainly by changes in stirring speed. 
An average of 0.55 w. per liter of medium was required by the stirrer. 


a-Amylases Production by Aspergillus niger PRL 558 

A medium containing 1.0% yeast extract, 2.0% starch, 0.5% calcium 
carbonate, salts (11), and 5% inoculum was fermented in the same manner as 
that for LU’. zeae except that the shaker speed was set at 260 r.p.m. and the 
fermentors were coupled at the beginning of the fermentation. 

The results are summarized in Fig. 12. The longer initial lag phase in the 
fermentor caused a considerable amount of overshooting of the stirring speed. 
However, a steady state was reached later. During this period the stirring 
speed varied between 310 and 340 r.p.m. and the rate of oxygen uptake was 
200 ml. per liter per hour. The average power required to maintain such an 
aeration was 0.4 w. per liter of medium. The a-amylase and cell yields obtained 
from the tank were comparable to those obtained from the shake flask. 

These results show that the oxygen uptake in a shake flask may be duplicated 
in a stirred tank with the help of the present translation device, and the 
corresponding aeration conditions for the stirred tank may be determined. 
One should be reminded that the master fermentor may be of any type, but 
the slave fermentor has to be of agitated type and its magnitude of aeration 
may be delicately controlled. Because of the limiting nature of the oxygen 
supply system, this device is only suitable for smaller fermentors with which 
extensive aeration studies are made. 
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THE STORAGE AND TRANSFER OF LOW POTENTIAL HEAT 
I. GLAUBER’S SALT AS A HEAT STORAGE MATERIAL! 


By J. W. HopGins anp T. W. HorrMan 


ABSTRACT 


A study has been made of the variables affecting the storage of low potential 
heat as the heat of fusion of Glauber’s salt. The geometry of the container was 
determined as the result of preliminary experiments, using a moving transfer 
surface. Extensive experiments with the final design yielded the following in- 
formation: 

(a) Low over-all efficiency is the result of segregation of phases as crystal 
growth proceeds. 

(b) Short heat paths in the exchanger are essential to keep thermal resistance 
at a minimum. 

(c) The crystal growth rate in this system is the rate governing process once 
phase isolation has occurred. 

(d) Over-all heat transfer coefficients are in the order of 2 B.t.u./hr. (sq. ft.) 
Pret 


INTRODUCTION 


The utilization of solar energy for domestic heating resolves itself into three 
separate problems: the collection of heat from the sun, the storage of the heat 
until required, and the employment of the stored heat for keeping a dwelling 
at a suitable temperature. 

The first of these problems has been exhaustively investigated by Hottel 
and others at the Massachusetts Institute of Technology (1) and Lof and 
associates at the University of Denver (2). The collector panel design, which 
represents a compromise between ideal performance and reasonable economy, 
in general consists of two or more glass plates separated by an air gap, with a 
black metal absorbing surface receiving the energy transmitted through the 
glass, which also reduces back radiation. The heat is removed from the black 
plate by a transfer fluid, usually air, and absorbed in the storage unit before 
the air is recycled to the panel. The optimum angle of tilt for collection of 
winter solar energy at the latitude of Kingston (44° 14’ north) is 62° from the 
horizontal facing due south, which will provide for 12°% better operation than 
will a vertical panel (3). 

The storage of the heat collected by the panel is, of course, of supreme im- 
portance to any solar heating system. The present investigation is an attempt 
to elucidate some of the controlling variables involved in the storage of low 
potential heat, and its transfer on demand. 

Although the transfer of heat from a panel to a storage unit will be efficient 
if the collector temperature is kept high, the efficiency of the acceptance of 
solar energy will diminish with rising plate temperature. Thus, if the efficiency 

‘Manuscript received April 1, 1956. 
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of the collection step is to be kept to a reasonable figure, it has been found 
necessary to maintain the collector plate below 120° F., which enables col- 
lection to be made of early morning and late afternoon radiation. In order to 
keep construction costs within reason, and because of low night temperatures, 
air has been chosen as the heat transfer medium, despite its low transfer co- 
efficients. The storage unit must accept heat from the air at 120° F. and must 
transfer it to a dwelling area which is about 70° F. 

Heat storage capacity can be provided in two ways—as heat capacity or as 
heat of fusion. Heat of fusion storage has already been used in this con- 
nection (5), but no detailed study has been made which will permit of the best 
design of a storage-transfer unit. The principal advantage of storing heat by 
chemicals is the saving in space—for example Glauber’s salt (Na2SO, - 10H,O) 
will store nearly five times as much heat as an equal volume of water between 
80-—120° F. Ice itself is the best storage material because it has the highest heat 
of fusion per unit weight or volume. However, because the melting point is 
not within the comfort range, a higher melting material must be used. The 
limitations on the storage material are fairly severe; it must be cheap to be 
practical, and it must have a melting and freezing point between 100° F. and 
85° F. The substances which best fulfill these criteria are crystalline hydrates, 
whose heat of fusion is higher, the greater the proportion of water of hydration 
in the molecule. 

All of the salts considered as heat storage materials undercool to a marked 
extent, and therefore nucleating material must be found which has the fol- 
lowing requisites. 

(a) It will act as a seeding material—that is, it must be nearly isomorphic 
with the hydrate used. 

(6) It must be only sparingly soluble in the solution resulting from the 
melting of the crystal. Any significant solubility of course lowers the melting 
point. 

(c) It must survive its environment through repeated cycles of melting 
and crystallization. 

In a heat storage apparatus, crystallization will start on the heat transfer 
surface where the nucleating material is present—that is, on the bottom 
surface in a stationary system where the nucleating agent is denser than the 
hydrate. The resulting crystal mantle will then present resistance to the flow 
of heat and to convection in the exchanger. In this investigation a study has 
been made of the effect of container geometry, agitation, and rate of heat 
withdrawal upon the operation of heat storage equipment with a view to 
providing design criteria for a full-size unit. All of the experiments were done 
using Glauber’s salt as the storage material, since a satisfactory nucleating 
agent was already at hand (6). It melts at about 90° F. with a heat of fusion 
of 104 B.t.u. per pound, and can be purchased for about $9.00 per ton. 

Dr. Maria Telkes has found that borax (Na2B,O; - 10H,O) is a satisfactory 
nucleating agent to induce crystallization of Glauber’s salt (6). At 90° F. the 
solubility is 2.71 gm. per 100 gm. of saturated sodium sulphate solution (4), 
at which concentration the phase transition temperature is lowered by about 
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Fic. 1. Phase equilibrium diagram for the system sodium sulphate — water. 
Fic. 2. Time-temperature cooling curves. 


2.5° F. In order to discuss the experimental results, it is appropriate at this 
point to examine the phase equilibrium diagram for the system sodium sul- 
phate — water (Fig. 1), and to consider the effects of cooling a system with 
concentrations of A, B, and C respectively. It is important to keep in mind that 
these considerations apply only to equilibrium conditions. 

Point A in Fig. 1 represents a system having the stoichiometric composition 
of Glauber’s salt, but at a temperature above the melting point of the hydrate. 
Because of the incongruent melting point, there is a considerable proportion 
of solid anhydrous Na2SO, present—a proportion which increases with rising 
temperature. Because of its high density, it rests on the bottom of the con- 
tainer. If one cools the mixture at A, constantly maintaining equilibrium con- 
ditions, the hydrated salt begins to crystallize at 90.3° F., and the temperature 
remains constant until the entire melt has gone solid, with the composition 
NaeSO,- 10H,O. The time-temperature cooling curve is represented schemat- 
ically in Fig. 2 (a). If the same cooling operation be performed on a mixture 
of composition B, the temperature halt at the phase transition temperature 
will be shorter, and at equilibrium, the two phases will consist of Glauber’s 
salt and a solution saturated with respect to Glauber’s salt, in the relative 
amounts e: f (Fig. 1). Further heat liberation will be possible only on lowering 
the temperature below 90.3° F. The cooling curve for this mixture is illustrated 
in Fig. 2 (6). Cooling of the saturated solution (composition A, Fig. 1), will 
not give any plateau at 90.3° F., and heat liberation will result only in con- 
tinual lowering of the temperature. Its typical cooling curve is illustrated in 
Fig. 2 (c). Therefore, maximum heat at the highest possible steady temperature 
is available when sodium sulphate and water are present in the stoichiometric 
ratio of 1: 10. 
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EXPERIMENTAL PROCEDURE 


The experiments were done in two stages. The first apparatus was designed 
to illustrate the efficiency of moving the heat transfer surface with respect to 
the crystal mass, in order to shed the crystal mantle which forms. This appa- 
ratus (illustrated in Fig. 3), consisted of a rotating coil in the crystal mass, 
(Glauber’s salt and 4 wt.-% of Na,B,O;- 10H2O) carrying water which re- 





Fic. 3. Rotating coil storage-transfer unit. 


moved the heat. An insulating Lucite box surrounding the jar containing the 
crystals was kept at the same temperature as the salt, thus permitting obser- 
vation without heat transfer through the container walls. The Glauber’s salt 
was first melted by hot water passed through the coil, and when the temperature 
of the mixture reached a suitable point, heat withdrawal was begun by cool 
water passed through the rotating coil. Temperatures of the water stream and 
of the salt mass were read on copper-constantan thermocouples. 

The second apparatus, built as a result of the experience with the first, had 
an entirely different form. In this exchanger, the salt was contained in a 
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shallow copper tray, 24 X 6 X 1 in., with air ducts of the same dimensions 
above and below the tray. With the exception of the tray itself, all of the 
components of this apparatus were thermal insulators—the outer walls of the 
ducts were made of asbestos board, and the entire apparatus encased in 
Foamglas insulation. The apparatus is illustrated in Fig. 4. Air, metered by a 
rotameter A, passed through a heating tube B into the exchanger. The first 
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Fic. 4. Shallow tray storage-transfer unit. 


four inches of the chamber were packed with sheets of glass wool air filter 
material to provide a calming section. Ahead of the copper tray a dummy 
section of insulating board (H) protruded six inches in front of the copper 
tray, to ensure stable boundary layer formation at the exchanger surfaces. 
Thus the tray section simulates the center portion of a practical heat exchanger. 
Air temperatures were measured by shielded copper-constantan thermo- 
couples, constantly monitored on automatic recording potentiometers. Since 
most of the air flowrates were in the streamline or transition region, it was 
necessary to use a three-pass baffle mixer just ahead of the air outlet thermo- 
couple. This caused impingement of the air on the copper sheet, and thus an 
increase in the transfer at that point. This spurious effect was minimized by a 
cemented asbestos sheet at the location of the impinging air. Altogether, 39 
thermocouples were used to measure air and plate temperatures and gradients 
in the solution. A calibrated hot wire anemometer was used to measure the 
air flow in the ducts, thus providing assurance that the mass flowrate of the 
air was the same in each duct. As in the first apparatus, the salt was first 
melted at an air temperature controlled by the thermoregulator, and heat was 
then withdrawn from the tray by cool air passed on either side. Here too, the 
filling was Glauber’s salt and 4 wt.-% of Na2B,O; - 10H.O. 
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RESULTS AND DISCUSSION 


In the first apparatus, crystallization started, as one would expect, upon the 
heat transfer surface of the rotating coil. Movement of the surface with respect 
to the saturated solution did nothing to remove the crystal mantle. Crystal- 
lization also occurred at the bottom of the container, and as the crystals grew, 
the hard crystal mass soon prevented coil rotation. The crystals in the bottom 
were obviously of two kinds—the hydrated crystals formed an impervious 
blanket over the anhydrous salt, thus stopping any subsequent hydration 
reaction. From this experiment the following conclusions were drawn: 


(i) Stratification of the nucleating agent does not seem to be a serious 
deterrent to crystallization. 

(ii) Crystallization occurs on the heat transfer surface, and the crystals 
are not detached by movement of the solution past this surface. 

(iii) The heat exchanger must be designed without internally moving parts. 

(iv) The hard crystal structure precludes the use of large tank-like con- 
tainers, because of the difficulty of transferring heat to melt the crystal mass. 

(v) The container must be shallow, so that the heat paths are short. 


All of the quantitative data were obtained using the second apparatus 
(Fig. 4). A typical experiment is illustrated by the time-temperature cooling 
curves shown in Fig. 5, in which an air flowrate of 12.6 c.f.m. (Re = 2300) was 
used. The flowrate was chosen as the result of considerations of the probable 
dimensions of a practical storage unit, and the heat requirement of 50,000 B.t.u. 
per hour—a normal value for a medium-sized dwelling in the locality of 
Kingston. The solution temperatures shown represent a portion of the solution 
which is unaffected by the impingement effect. The bottom plate temperatures 
at various locations are shown as the run proceeds. The undercooling amounted 
to 6.3° F., at which time crystallization began. There is no appreciable plateau 
at the maximum, and the solution temperature quickly reaches a lower limit 
of 81° to 82° F. Air at lower temperatures than this will be useless for space 
heating. 

Heat was automatically turned on when the air temperatures reached 82° F., 
and the heating of the tray began. A short temperature plateau was observed, 
followed by a sharp rise; the melting point is markedly higher than the highest 
temperature reached during heat extraction. At any particular time, the 
quantity of heat transferred is 


dq = W. Cp. dd. (Tout — Tin), where 

Cp = specific heat of the air, B.t.u./Ib. (°F.). 
dé = time interval, hr. 

W = flowrate, in lb./hr. 


By integrating under the time-temperature curves, the total heat trans- 
ferred over a given period was calculated; from these measurements, the total 
heat-in was found to equal the total heat-out, and this heat balance indicates 
negligible heat loss in the exchanger. Over the crystallization period, 94 B.t.u. 
of heat were realized, as compared with the total theoretical amount of heat 
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Fic. 5. Time-temperature cooling curves for tray exchanger containing Glauber’s salt. 


capable of being stored in the test section, 570 B.t.u. This same efficiency of 
16.5% was obtained over the range of air flows used—from Re 1000 to Re 2300. 

In Fig. 6 is shown a piece of crystal taken from the exchanger after the runs 
were complete, which illustrates the two phases, anhydrous sodium sulphate, 
surmounted by Glauber’s salt. The hydrated salt so effectively isolates the 
anhydrous layer from the solution that the supernatant liquid never crystal- 
lizes. Because of the inverse solubility of the anhydrous salt (it becomes less 








300 CANADIAN JOURNAL OF TECHNOLOGY. VOL. 33 


soluble with rising temperature), the crystals form a very dense cake on 
repeated cycling of temperature. That is, on heating, additional anhydrous 
material deposits in all the crevices, while on cooling these deposits are only 
slightly dissolved before the isolating hydrate mantle forms. 





Fic. 6. Crystal formation from tray heat exchanger. 


Once the blanket of Na2SOQ,-10H,O has isolated the anhydrous salt, the 
system acts as a saturated solution, which explains the absence of a long 
plateau in the cooling curve. Further crystallization will occur only as a result 
of lowering the temperature. Since rate of crystallization is a function of the 
degree of undercooling, equilibrium conditions do not prevail once this isolation 
is complete, and the temperature attained is then a function of the rate of 
withdrawal and the inherent crystal growth rate. 

The 16.5% efficiency observed in these experiments represents the worst 
efficiency. In practical use the withdrawal of heat would be intermittent, and 
during the off cycle the equilibrium conditions would be at least partially 
regained. Nonetheless, it would still behave as a saturated solution, for no 
reaction with the anhydrous salt on ‘the bottom is possible. 
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Fig. 7 illustrates the variation with time of the total heat transferred and 
the over-all heat transfer coefficient. The large heat transfer at the beginning 
of the cooling and heating runs represents the heat capacity of the apparatus. 
The minimum in the cooling run denotes the beginning of heat generation 
and the decreasing of internal thermal gradients. As the plate temperature 
increases, the over-all heat transferred from the bottom plate reaches a maxi- 
mum, while the over-all heat transfer coefficients remain fairly constant, 
decreasing only slightly as the salt builds up on the bottom. Since there is no 
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Fic. 7. Heat transfer coefficients and heat transferred in tray heat exchanger. 
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crystal growth on the top plate, heat reaches it only by convection and con- 
duction through the liquid, resulting in a lower heat transfer. Because the 
nucleating agent is denser than its environment, crystallization starts at the 
at the top of the insulating layer of anhydrous sodium 





worst possible place 
sulphate. 

On the heating cycle, the total heat transferred and the heat transfer co- 
efficient for the bottom plate reach a sharp maximum. This maximum repre- 
sents the point where the internal surface contact changes from solid to liquid, 
allowing for better heat acceptance. Heat transferred through the top plate 
is low, for it is removed from the plate by conduction through the liquid. The 
plate reaches a high temperature quickly, and therefore high internal tempera- 
ture gradients are set up. 


CONCLUSION 

Although Glauber’s salt has been used in an actual heating installation (5), 
the present investigation has indicated that it is not a very satisfactory storage 
material. The manner in which equilibrium is hindered by the growth of the 
hydrate to isolate the anhydrous salt lowers the efficiency at the outset; and 
the fact that crystal growth from the saturated solution cannot keep up with 
the rate of heat abstraction means that the slow step in the heat removal cycle 
is governed by the inherent crystal growth rate. Other more promising systems 
are under current investigation. The features being sought are homogeneity 
above the melting point, ease and reproducibility of nucleation, and rapidity 
of crystal growth. 
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AN OPERATION-TIME CONTROLLER FOR A RECORDING POTENTIOMETER' 


By G. O. HANDEGORD AND D. G. COLE 


Certain applications of recording potentiometers require that a record be 
made for only part of each day. In some cases, such as solar radiation measure- 
ments, the ‘‘on’’ period of the recorder does not coincide with normal working 
hours and manual starting and stopping of the instrument is not possible. 
The instrument can be operated continuously but this results in an excessive 
amount of chart being used, particularly in the winter months when sunshine 
occurs for only eight hours. 

A standard time switch, such as the ‘‘Inter-matic’” Model TS-60SP, manu- 
factured by the International Register Co., may be employed to stop and 
start the instrument at any predetermined time. This time switch is commonly 
used in business places to operate display lighting during the dark hours of 
the evening. 

When the potentiometer used incorporates a standard cell circuit with 
automatic standardization, one difficulty arises in that the time switch may 
stop the instrument when the standardizing circuit is closed, resulting in an 
undesirable drain on the standard cell during the ‘‘off’’ period. This difficulty 
may be overcome by incorporating a micro switch in parallel with the time 
switch to operate the instrument should the time switch open while standard- 
ization is in progress. 

The Brown 153XI1IV-X-27, strip chart recording potentiometer used for 
solar radiation measurements in Saskatoon has been modified in the above 
manner. The time switch and micro switch are incorporated into the instru- 
ment as shown in Fig. 1. The micro switch employed is a normally open type 
#BZ-2R, installed so that its contact pin bears against the standardization 
slide wire mechanism. During standardization, the mechanism is moved in 
towards the chassis so that the micro switch pin protrudes fully, thus closing 
the switch and operating the instrument until standardization is complete. 
The micro switch is mounted on a brass plate which is secured to one of the 
mounting posts of the standardization assembly cover as shown in Fig. 2. 
Both the time switch and micro switch are wired in parallel with the drive 
switch as shown in Fig. 3. This wiring arrangement does not involve any change 
in the original wiring of the instrument and all the connections may be made 
to the terminal board on the rear of the chassis. 

The time switch and micro switch operate only the chart drive circuit of 
the instrument, power being supplied continuously during both ‘‘on’’ and“‘oft"’ 
periods to the amplifier circuit as recommended by the manufacturer. Inter- 
mittent operation of the recorder as determined by the setting of the time 
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Fic. 1. Location of time switch and micro switch in recorder case. 
Fic. 2. Close-up view of micro switch and mounting bracket. 
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switch is achieved with the instrument chart drive switch in the ‘‘off’’ position. 
Continuous operation of the recorder is readily obtained by moving the chart 
drive switch to the ‘‘on”’ position. 
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Fic. 3. Chart motor wiring diagram. 


The total cost of materials required for this automatic switching arrange- 
ment is approximately thirteen dollars. The savings involved in chart con- 
sumption for the particular application cited will be approximately eighteen 
dollars per year, based on an average period of 14 hr. of sunshine and a chart 
cost of four dollars per roll. 
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